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, Gravity (QEG) on an infinite dimensional theory space. We consider "conformally 

reduced" gravity where only fluctuations of the conformal factor are quantized and 
^ employ the Local Potential Approximation for its effective average action. The 

' requirement of "background independence" in quantum gravity entails a partial dif- 



ferential equation governing the scale dependence of the potential for the conformal 
factor which differs significantly from that of a scalar matter field. In the infinite di- 
mensional space of potential functions we find a Gaussian as well as a non-Gaussian 
fixed point which provides further evidence for the viability of the asymptotic safety 
scenario. The analog of the invariant cubic in the curvature which spoils pertur- 
bative renormalizability is seen to be unproblematic for the asymptotic safety of 
the conformally reduced theory. The scaling fields and dimensions of both fixed 
points are obtained explicitly and possible implications for the predictivity of the 
theory are discussed. Spacetime manifolds with R'^ as well as S'^ topology are con- 
sidered. Solving the flow equation for the potential numerically we obtain examples 
of renormalization group trajectories inside the ultraviolet critical surface of the 
non-Gaussian fixed point. The quantum theories based upon some of them show a 
phase transition from the familiar (low energy) phase of gravity with spontaneously 
broken diffeomorphism invariance to a new phase of unbroken diffeomorphism in- 
variance; the latter phase is characterized by a vanishing expectation value of the 
metric. 



1 Introduction 



The formulation of a mathematically consistent quantum theory of gravity whose 
predictive power extends down to arbitrarily small distances continues to be one of the 
major challenges of modern theoretical physics [1]. Even though the first attempts at 
quantizing the gravitational field date back to the 1930s [2] the various approaches avail- 
able today are still plagued by considerable technical and conceptual problems. On the 
conceptual side, the "background independence" [3-5] of classical General Relativity and 
of the yet to be constructed theory of quantum gravity is at the heart of many difficulties. 
In the asymptotic safety approach to quantum gravity [6-29] the gravitational degrees of 
freedom are described by a metric Qf^^, as in General Relativity, and the standard rules of 
quantum mechanics are applied to it. Nevertheless, because of the requirement of "back- 
ground independence", the resulting quantum field theory of gravity has many features 
which are rather different from the familiar matter field theories on a fixed, non-dynamical 
spacetime. These unusual features emerge because the theory is supposed to dynamically 
explain the origin of the spacetime it "lives" on as a kind of ground state property, and 
as a result, in its construction no metric should be distinguished from any other a priori. 
Most investigations of the asymptotic safety scenario are based upon the effective average 
action for gravity [7] . In this approach the requirement of "background independence" in 
the above sense of the word is met, somewhat paradoxically, by means of the background 
field technique [30]. Even though, formally, the gravitational average action obtains by 
a similar Wilsonian mode cutoff and coarse graining operation as its counterpart for or- 
dinary matter field theories [31-34] there is a conceptually crucial difference: The metric 
determines all physical (proper) length and mass scales; in particular it fixes the phys- 
ical coarse graining scale which is to be attributed to the cutoff parameter k built into 
the average action Ffe. Hence in quantum gravity the problem is that the field that is 
"averaged" or "coarse grained" by itself has to define the physical meaning of this coarse 
graining operation. If one was to use a rigid metric for this purpose the parameter k would 
have no physical interpretation; hence Fj. hardly could be used for "renormalization group 
improving" classical solutions or field equations [35-47], for instance. 
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independent"!^ 



In a recent paper [22], henceforth referred to as [I], we explained in detail how the 
use of the background field technique, besides making a diffeomorphism invariant func- 
ike the average action and its renormalization group (RG) flow "background 
in the sense that nowhere in the construction any preferred (rigid) metric 
occurs. In order to disentangle this problem from the gauge fixing issues which are the 
usual motivation for the background method we considered the "conformal reduction" 
of Quantum Einstein Gravity (QEG) in which only the conformal factor of the metric 
is quantized and all other degrees of freedom are discarded. Performing the conformal 
reduction on top of the usual Einstein-Hilbert truncation one is led to quantize an ac- 
tion of the familiar 0^-type. As we saw in [I] the RG flow in this conformally reduced 
Einstein-Hilbert (or "CREH" ) truncation is different from that in a standard scalar mat- 
ter field theory, the reason being that in the gravitational case the mode suppression 
term exp(— A^S") in the functional integral defining must have a specific background 
field dependence in order to comply with the requirement of "background independence" . 
This extra field dependence of AkS is absent in the matter field theory. As a consequence, 
gravity in the CREH approximation has an RG flow which is very different from that of 
a standard 0^-theory. Including the extra field dependence turns the merely logarithmic 
running of the cosmological constant (earlier discussed by Polyakov [51] and by Jackiw et 
al. [52]) into the expected A;*^-running, for instance. (Similar observations had been made 
by Floreanini and Percacci [53] in a different theory.) 

Remarkably, the RG flow of the CREH approximation turned out to be qualitatively 
identical to the full Einstein-Hilbert flow where all degrees of freedom are quantized. In 
particular it has a non-Gaussian flxed point (NGFP) suitable for the asymptotic safety 
construction. The counterpart of the CREH approximation on the side of the standard 
matter theories is Symanzik's asymptotically free 0^-theory with a negative quartic cou- 
pling [54,55]. This theory has a potential which is unbounded below, and it is notorious 



^Here and in the following the term "background independent", put in quotation marks, means the 
absence of a preferred metric; in this sense the term is frequently used in loop quantum gravity [3-5] or in 
the dynamical triangulation approach [48-50] , for instance. Referring to the background field formalism, 
no quotation marks will be put. 
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for the resulting infrared instability. In the context of QEG this instability is not an issue, 
however. In the full theory it is presumably cured by higher derivative terms. In fact, it 
is known [11] that a J y/g' R'^-term, added to the Einstein-Hilbert action, stabilizes the 
kinetic term of the conformal factor at high scales; cf. also ref. [56]. 

The upshot of the analysis in [ I ] was that the quantum fluctuations of the conformal 
factor alone seem to be typical of the full set of degrees of freedom in g^^, as far as their 
impact on the RG flow is concerned. This suggests that also in more general truncations 
the conformal factor should play a representative rolqj. Since the scalar-like reduced 
theory is technically much simpler than full-fledged QEG it suggests itself as a tool which 
should allow us to probe regions of theory space which are too difficult computationally 
in the full theory. 

In the present paper we start this program with the local potential approximation, 
or "LPA", which includes arbitrary non-derivative terms in the truncation ansatz. In 
fact, the CREH truncation (in 4 Euclidean dimensions) obtains by inserting metrics of 
the form (yf^^, = 0^(x) ^^i, into the Einstein-Hilbert action. Here ^f^^ is a fixed reference 
metric. The result is a kinetic term for which, apart from its sign, is of the standard 
form oc d^(j) d^(f), as well as a potential with two terms, proportional to 0^ and 0^, 
respectively: 



CREH 



k 



(0) = -^ (3^i?0^-|A,0^). (1.1) 



(Here R is the curvature scalar of ^^,y.) The idea of the LPA [33] is to retain the standard 
kinetic term from the CREH approximation but to allow for arbitrary non-derivative 
terms: 



Tk[(p = const] = + Jd^x^/f ^7^(0). (1.2) 



The running potential Uk{4>) may develop any 0-dependence; it is no longer restricted to 
be a quartic polynomial. The LPA has played an important role in the RG analysis of 
standard scalar theories, in particular in the context of spontaneous symmetry breaking 
and the approach to convexity [33]. It is therefore the obvious next step in the analysis 



^For earlier work in a similar spirit see ref. [57]. 
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of conformally reduced QEG. The relevant theory space is infinite dimensional now; the 
corresponding fiow equations will include a partial differential equation for Uk{ ■ )■ 

Clearly the LPA cannot be a numerically precise approximation to the full theory, 
even if the fiuctuations of have always a "typical" impact on the RG fiow. Nevertheless 
it should provide us with some qualitative understanding and conceptual insights into the 
general properties of QEG. 

The (relative!) technical simplicity of the reduced theory comes at a certain price, 
however. In general it will not be possible to relate terms in Uk{(f)) to invariants I[giiu] 
depending on the full metric in a unique way. Typically there will be many invariants 
contributing to a given term in the potential when we insert g^i, = cfP' ^^y. Nevertheless, 
the LPA is sensitive to a certain projection (or "shadow") of the invariants I[g^u] which 
make up the general Tklg^^^]. If g^^, is the metric on a curved space, Uk{(p) encodes 
information about higher derivative terms, for instance. If we consider a series expansion 
in the curvature scalar, for instance, 

^k[g,.] = J2crXk) fd'x^Rigr (1.3) 

and choose g^y to be the metric on a unit 4-sphere, then the corresponding potential 

Tfc [0 = const] I J d^x reads 

f/fc(0) = ^(12)"c„,(A;)0^-2". (1.4) 

n 

The constant term in (11.41) . n = 2, descends from the invariant J d^x^/g' R^, but also 
/ d'^Xy/g' {Rf^uY or f d^x^/g {R^^apY would contribute to the constant piece in Uk{4>), as 
would many non-local terms. Similarly, the monomial receives contributions from R^ 
and the Goroff-Sagnotti term J d'^x y/g' Rp.uap R""^ R'"^^^ [58], among others. 

Despite these limitations the experience with conventional field theories suggests 
that the LPA is a powerful tool for a first exploration of regions in theory space which are 
computationally inaccessible otherwise. In the present paper we shall employ it in order 
to analyze and illustrate several conceptual issues related to the renormalization group 
fiow of Fjt which are as to yet too computationally demanding to be studied in the full 
theory. In particular we focus on the following three topics. 
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(a) In the LPA we have the possibihty to obtain information about the non-Gaussian 
fixed point of the gravitational average action on an infinite dimensional theory space. 
We shall find that the RG fiow on the space of (dimensionless) potential functions and 
Newton's constant admits both a Gaussian and non-Gaussian fixed point. We compute 
the fixed point potentials and perform a linear stability analysis in each case. The scaling 
dimensions of the Gaussian fixed point (GFP) will be seen to characterize the difference 
between a standard scalar and the conformal factor in a particularly clearcut way. The 
scahng fields and dimensions (critical exponents 9) of the NGFP are of interest for the 
asymptotic safety construction. 

The relevant scahng fields, those with Re^ > 0, decay for A; — > oo, while the 
irrelevant ones, with Re^ < 0, blow up in this limit. An asymptotically safe quantum 
theory is specified by any RG trajectory that hits the NGFP for k oo. We define 
the ultraviolet (UV) critical manifold ^uv to be the set of all points in theory space 
which arc pulled into the NGFP under the inverse RG flow (i.e., for increasing k). In 
the simplest case ^uv is a smooth submanifold of theory space. Every complete RG 
trajectory which lies inside =5ljv defines a possible quantum theory. Hence, with Auv = 
dim =5ijv the dimensionality of the critical manifold, there exists a Auy-parameter family 
of asymptotically safe quantum field theories; which one is actually realized in Nature 
must be decided by an experimental determination of exactly those Auy-parameters whose 
values are not predicted by theory. If t5^uv is indeed a smooth manifold, and the NGFP has 
no marginal scaling fields (Re^ = 0), then Ayy equals the dimensionality of a subspace 
of the tangent space to ^uv at the NGFP, namely the one spanned by the relevant 
scaling fields, i. e. those whose critical exponents 9 have a positive real part. Therefore 
the number of relevant scaling fields is related to the degree of predictivity that can be 
achieved. (See [29] for a more general discussion.) 

Within the LPA we will be able to study ^uv near the NGFP in great detail. We 
shall find that the value of Auv depends sensitively on the space of functions in which 
allowed scaling fields are supposed to "live" . It turns out that in principle there could be 
infinitely many relevant ones, but as we shall see this question cannot be decided on the 
basis of the LPA alone. 
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(b) Using numerical methods we shall analyze the RG evolution of Uk{4>) also in the 
nonlinear regime. We are particularly interested in phase transitions, triggered by a 
variation of k, in which the global minimum of either jumps, or continuously evolves 
from a value 7^ to = 0. We interpret this as a transition from the familiar phase of 
gravity with broken diffeomorphism invariance to a phase with unbroken dijfeomorphism 
invariance. 

We are considering metrics of the form g^^^, — cf)^ 'g^i, where (j) is the expectation 
value of the quantum conformal factor and g^i, is a non-degenerate classical reference 
metric. If the expectation value 0, interpreted here as the minimum of C/^, is nonzero the 
expectation value of the metric, g^i, = (7^;/) , is non-degenerate as well. But if = then 
the metric has a vanishing expectation value. (More precisely, in view of g^i, — Sab 
and since g^i, = 0^ depends on (p quadratically, it is the vielbein that has a vanishing 
expectation value.) 

The situation where = can be thought of as the restoration of a spontaneously 
broken symmetry [53,59]. The relevant local symmetries are the spacetime diffeomor- 
phisms. If the quantum state is such that {'j/xi^) 7^ the group of diffeomorphisms is 
spontaneously broken down to the stability group of (7^1/). If (7^1/) = ?7^i/, say, this sta- 
bility group is the Poincare group. The broken phase of gravity is the one we are familiar 
with. It is conceivable though that the symmetry breaking effectively disappears at very 
small distances. This is exactly what happens in spontaneously broken Yang-Mills theo- 
ries, in the standard model, say. If we study phenomena involving momentum scales far 
above the W-mass, the full SU{2) x U{1) symmetry appears unbroken. Technically this 
symmetry restoration can be inferred from the running Higgs potential uf^^^^ at k ^ M^. 

We are going to analyze the running potential of the conformal factor in the same 
spirit. We shall find that typical trajectories in S^uy, for k sufficiently large, describe a 
phase of unbroken diffeomorphism invariance in the sense that the metric has a vanishing 
expectation value there. 

(c) While the beta functions implied by the functional RG equations do not refer to any 
special metric they do depend on the topology of the manifold which carries the metrics 
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g^ivi g^yi and 'g^^. We illustrate this phenomenon by discussing in parallel the two 
cases where the topology of the (Euclidean) spacetime manifold is that of a sphere 5*^ 
and of flat space i?^, respectively. The beta functions (or rather, functionals) of Uk{-) 
are different in the two cases since they depend on the spectrum of the Laplace-Beltrami 
operator built from the corresponding ^^i,. 

The remaining sections of this paper are organized as follows. In Section [2] we 
review and expand on conformally reduced QEG; in particular we show how the back- 
ground field method achieves "background independence", and we introduce two types 
of background-quantum field split symmetries. The LPA is introduced in Section [3] then 
and the corresponding flow equations are derived for the 5*^ and the topologies. In 
Section m we analyze various general properties of the LPA flow, and Section [5] is devoted 
to its fixed points; for both the Gaussian and the non-Gaussian fixed point we determine 
the fixed point potential and perform a linear stability analysis in order to determine the 
local tangent space to the UV critical manifold. In Section [6] we solve the partial differen- 
tial equation for the running potential numerically and describe various scenarios for the 
phase transition to a phase of (conformally reduced) QEG with unbroken diffeomorphism 
invariance. The Conclusions are given in Section [71 



2 "Background Independence" via Background Fields 

In the construction of the full-fledged gravitational average action [7] there are three 
different metrics that play a role: the microscopic metric, i. e. the integration variable in 
the path integral, 'j^^, its expectation value g^i, = (7^^), and the background metric 
In conformally reduced gravity [I] all three of them are assumed conformal to the same, 
fixed reference metric 'g^^. We write the conformal factor of 7^1, and g^^, as x^'^i^) and 
Xb'(^) where, in d dimensions, u = 2/{d — 2). The motivation for the exponent 21/ is that 
in this way the Einstein-Hilbert action provides a standard kinetic term for x- 

In the approach of [ I ] the original path integral over 7^,^ is reduced to an integral over 
x{x) only, and the resulting scalar-like theory is then reformulated using the background 
field method. The "microscopic" conformal factor x is decomposed as x = Xb + / and 
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the x-integral is replaced by an integral over /. We denote the expectation values of / 
and X by 7 = (/) and 

0=(x) = XB + 7, (2.1) 
respectively. Thus the foUowing metrics are to be distinguished: 

= X^"" g,u = (XB + ff'' 9,. (2.2) 

V ^ di^i^ (2.3) 

9,u = {iXB + ff') 9,u (2.4) 

gtMi^ = ^^'^ Qnu = (xb + {DY'^ g„v (2.5) 

Since 2z/ 7^ 1 in general, wc arc performing a nonlinear background split here which causes 
the fields g^j^u and gjj,^ to be different. In this respect the correspondence between the full 
and the conformally reduced theory is not perfect. In the full theory one performs the 
linear background split — g^^^ + h^^, integrates over h^,,, and defines /i^,^ = (h^,^). In 
the latter case there is no distinction between g/j^i, and ^^j,: 

g^u = = g^,^ + (V) = ^M^^- (2-6) 

In [I] we explained in detail why, both in the full and the reduced theory, the 
background field is indispensable if one wants to implement "background independence" 
in the sense of [3-5] and give a physical interpretation to the RG parameter k. One 
expands the path integral in eigenmodes of the covariant Laplacian □ which is constructed 
from g^j^ and then introduces A; as a cutoff in the spectrum of □. In this manner the 
typical length scale displayed by the D-eigenmode with eigenvalue — A;^, the "last" mode 
integrated out, is an approximate measure for the proper (w. r. t. ^^^) coarse graining 
scale. The infrared (IR) cutoff is built into the path integral by adding an appropriate 
cutoff action A^S*!/; Xb] to the bare action. It is a quadratic form in / which involves a 
XB-dependent integral kernel 7^fc[xB]. (See [I] for further details.) 

In the full theory, for vanishing ghost fields, the gravitational average action is a 
functional of the expectation value and the background field, respectively: Tj^f/i^,^;^ ] = 
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^k[gfiu,gij,u]- III the reduced setting, depends likewise on / and xb or alternatively on 
(f) = Xb + f and xb- ^k[f]XB] = Tk[4>:XB]- Even though it is kept fixed usually, it is 
sometimes important to keep in mind that also depends on the reference metric g^,y. So 
occasionally it will be helpful to include it in the list of arguments, writing Tk[f; XB,dfMu] 
or Tk[(p,XB,g^lu], respectively. 

In [ I ] we derived the following FRGE for the reduced average action: 



kdkTk[f;XB,gf,u] = ^Tr 



^k\f'^XB,gf,u] +'R■k[xB,g^^u]) kdkTlkixB.g^tu] 



(2.7) 



This equation has the same structure as in a scalar theory, except that the properties of 
the cutoff operator TZ^ are different here. In [I] we explained in detail how it has to be 
constructed in order to be "background independent" and to take account of the fact that 
the metric itself determines the proper coarse graining scale. 

The reference metric 'g^v is unphysical and, in fact, has no counterpart in the full 
theory. A Weyl rescaling of can always be absorbed by a suitable redefinition of the 
conformal factor. This is formalized by saying that the "physical" metrics 7^1,, ^^j,, g^^, 
and g^y are invariant under the following background-quantum field split symmetries: 

g',u = e-'"^^^) g,u (2.8a) 

/' = e-^^^/^ /, Xb = e'^^^^/'^ Xb (2.8b) 

J' = g-(x)/.^ 0/ _ ^a(x)/^ ^_ ^2.8c) 

Since the decomposition of -y^i, as a conformal factor times 'g^^, is completely arbitrary, 
the quantum theory of the conformal factor must respect the split symmetry (12. 8p . In 
particular we require the cutoff action AkS[f;xB,^iJ.u] to be invariant under the "cr- 
transformations" (12. 8p . as we shall call them. This invariance has an important con- 
sequence for the solutions Tk[f; Xb, dfiu] of the FRGE (12. 7p : If F^ is invariant at one value 
of k, it is so also at any other k. The flow does not generate non-invariant terms and, as 
a consequence, it is enough to include cr-invariant terms in the truncation ansatz. 
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Later on we shall fix ^f^^, to be the metric of a round 4-sphere with radius r, S'^(r ), 
so that the corresponding line element reads, using standard coordinates, 

Qf,^ dx'' dx" = r 

If we insist on this form of the reference metric, the invariance under local cr-transforma- 
tions is broken to a symmetry under global ones only. Global split transformations (12.81) 
with, by definition, a{x) = a = const, preserve the form of this metric, but rescale the 
radius of the sphere according to 

r' = e-" r. (2.10) 

Once the conformal reduction is performed, we are left with a theory of the field 
x{x)- To quantize it, we perform a second, this time linear background-quantum field split 
by decomposing X = Xb + f ■ This decomposition is supposed to mimic the = 'g^y + h^j^y 
decomposition in the full theory. Obviously x ^"^^ = Xb + / are trivially invariant under 
a second set of split transformations: 

f = f + r{x), f = J+r{x) 

(2.11) 

Xb = Xb-t{x), 9!,^ = %u- 

While the "cx-transformations" have no analog in the full theory, the above "r-transfor- 
mations" do have a counterpart: h'^^ = h^^ + t^^, h'^^ = h^^ + r^^, g'^^ = g^^ - t^^. 

It is to be emphasized that both in the full and the reduced theory A^S* and as a 
consequence the FRGE and its exact solutions Tk are not invariant under the respective 
r-transformations if A; 7^ 0. This implies that the background field cannot be "shifted 
away" by redefining the quantum field (integration variable). In a certain sense, the 
background field must have acquired a "quasi physical" meaning therefore. From our 
discussion in [ I ] of the special status of the metric, it defines all proper length and mass 
scales, including that of it is clear what this meaning ijf|: The background encodes the 
information about the physical interpretation of the mode cutoff, the proper length scale of 

■^The gauge fixing issues we leave aside here. They are not important for the present discussion. 



dC^ + sin2 C 



drf + sin^ rj 



[dO 



2 I -2 

+ sm 



edct?) 



(2.9) 
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the "coarse grained" domains in spacetime, in particular. Being a cutoff in the spectrum 
of □, the parameter is a "proper" (as opposed to coordinate) length with respect 
to the background metric. The r-symmetry is broken by the cutoff action A^S which 
implements the mode suppression. For k ^ the functional approaches the standard 
effective action F. In this limit F^ becomes invariant under the "r-transformations" . 

In [I] we solved the FRGE (12.71) on the 2-dimensional truncated theory space 
spanned by the ansatz 

r47;xB,?^.] = -^^ Jd'x^/f hixB + 7) (-3 + 1^) (xb+7) 

(2.12) 

-|Afc (xB+7)' 

This functional is invariant under x-dependent a- as well as r-transformations. The 
pertinent cutoff action AkS[f;xB,^fiu] is only cr-invariant, however. A generalization of 
(I2.12P will be encountered in the next section. 

3 The Local Potential Approximation 

Next we introduce the local potential approximation (LPA) to the conformally re- 
duced theory of gravity discussed in the previous section and we derive the corresponding 
beta functions. We restrict the discussion to the case of a maximally symmetric reference 
metric ^^jy. To start with we shall assume that it refers to a sphere S'^(r ); later on we 
shall also consider the flat space R^. For clarity's sake we focus on d = 4 spacetime 
dimensions here; the d-dimensional generalizations of various key formulas are collected 
in the Appendix. 

3.1 The Truncation Ansatz 

The LPA truncation ansatz has the same kinetic term as the CREH truncation 
(I2.12P but instead of the single quartic term proportional to A^ we now allow for an 
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arbitrary running potential Hk'. 



))] 



3 



d'xVf \^UxB + f) (-3 + 1^) iXB + f) 
+ r''Hk{rixB+7))]- 



(3.1) 



After the restriction to spherical reference metrics we must make sure that Tk is invari- 
ant under x-independent a-transformations (12.81) which change the radius r according to 
f l2.10p . This explains the explicit factors of r in (13.11) . They guarantee that is indeed 
invariant under global cr-transformations: 



The ansatz (13.11) happens to be invariant under the r-transformations (12. lip . This 
accidental r-invariance is due to the still comparatively simple form of the truncation 
ansatz; recall that the FRGE and the exact solutions do not have this symmetry. In 
fact, in a more general truncation we would assume that the potential depends on the 
two fields Xb = ^Xb and / = r / separately, and not only via their sum (p = rcf). 
Then Tk would not be r-invariant in general. The simplifying assumption of a single- 
field LPA potential parallels analogous truncation steps in the full theory [7]: Starting 
from a general potential iffc(0, xb), without loss of generality, we may decompose it as 
Hk(4>,XB) = HkQ)) + HkQi.XB) where iJfc(0) = i^fe(0,0), and Hk vanishes for = xb- 
Adopting the ansatz (13.10 amounts to neglecting if^, i. e. to the approximation ^ Hk- 
With the potentials Hk and Hk corresponding to the functionals rfc[5f^,^] and Tk[gij.u,'gf^i^], 
respectively, this is the same kind of "single-field approximation" one often makes in 
full QEG, see Section 3 of ref. [7]. However, even though we take Hk to depend on one 
field only, special solutions of the type (13.10 are in fact general enough to demonstrate 
the qualitative consequences of the r-invariance breaking xe-dependences which enter the 
flow equation via TZk and are dictated by "background independence" , see [ I ] for a first 
example. 

It will be convenient to combine the i?-term which stems from J y/g^ R with the 
potential Hk- On the 4-sphere R = 12/r^ is a constant. Therefore, if we introduce the 



[e'^ /; e'^ XB,g,u {S\e-'' f ))] = /; xb,9,u {S\r 



))]• 



(3.2) 
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total potential as 

FS)=^' + Hk{'^) (3.3) 
the average action assumes the form 

+ f-^Ffc(f(xB+7))}. 



(3.4) 



We are going to insert this ansatz into the FRGE (12.71) . Under the trace on its RHS we 



need the Hessian \ It reads 



.(2) 
k 

^(2) FT. 



-□. + f-2F;(f(xB(x)+7(x))) 



1 xy 

3 



S{x - y) ^3.5) 



X] 



Here and in the following a prime denotes a derivative with respect to the argument. 
Comparing ([32]) to ([L2D we see that Uk{4>) = -^^^ f Ffc(f 0). 

The ansatz (13. 4p contains two scale dependent parameters, a single constant, namely 
Gk, and a function of the field, -Ffc( ■ )• We determine their respective beta functions by 
evaluating both sides of the FRGE with the ansatz inserted for special field configurations 
which project out the relevant monomials in the action. In order to find kdkFk{-) we 
set / = and xb = const, compute the trace on the RHS of the FRGE as a function of 
the x-independent variable xb, and equate the result to kdkFkixB)- The beta function 
for Newton's constant obtains from the kinetic term. We set xb = const and leave / 
nonzero and x-dependent. We expand the trace to second order in / and second order in 
the derivatives; the (9/)^-piece is proportional to kdk (l/Gk) then. 

As for both of these calculations a constant background field is sufficient we assume 
Xb{x) = const = Xb from now on. 

Before continuing a remark about dimensions might be in order here. In the above 
equations we use units such that the coordinates x'^ are dimensionless. Hence we have 
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the canonical mass dimensions 

[7m^] = [gf.u] = [^^J = [dt^u] = -2 

[/] = [7] = [Xb] = M = (3.6) 
[Hk] = [Fk] = -2, [r] = -1. 

3.2 The Cutoff Operator 

In [ I ] we explained in detail why it is crucial to define as a cutoff in the spectrum 
of □ rather than □, say. (Here □ and □ are the Laplace-Beltrami operators related to 
g^^ = Xb Qpiv S'lid g^u, respectively.) At least for truncations of the type considered here 
this can be achieved by designing the cutoff operator IZk in such a way that upon adding 
it to r^^"* is leads to the replacement 

(_n)_(_n) + fc2^(o)(_n). (3.7) 

Since, for xb = const, □ = Xb^ this substitution rule amounts to 

i-a)^{-n) + xlk'R^'\-^). (3.8) 

In our case T^^ is given by (13.51) . Therefore we define TZk such that, in operator notation. 
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As a result, we must set 



-□ + xl (-^) + F;(f (xB + / (x)))l . (3.9) 



n.lxB-, 9,.] = xl R^'^ (-^) . (3.10) 

Note the factors of Xb the RHS of eq. fl3.10p . They are typical of conformally reduced 
gravity and would not appear in the standard quantization of a scalar field. Besides the 
different transformation behavior under diffeomorphisms, they are the most important 
feature which distinguishes the conformal factor from scalar matter fields. They lead to 
very significant modifications of the RG fiow [ I ] . The cutoff action corresponding to ( ]3.10p 
reads: 

A,S[f;xB,9,.] = -^^ jd'x^ f{x)R^'\-£j^)f{x). (3.11) 
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This functional is indeed invariant under global cr-transformations, as it should be: 

AkS[e'^ f-e'^ XB,e-^'^g,,] = AkS[f;xB,g,u]. (3.12) 
Note, however, that AkS is not invariant under r-transformations: 

AkS[f + t;xb -r,g^,u] ^ 5* [/; xb, (3.13) 
At this point the flow equation assumes the form 
kdkTk[f;XB = const, gf,^{S^{r))] 



Tr 



-□ + xi (-^) + f FHinxB +7)) 



(3.14) 



y<G,kd,{G^'xlk'R^'\-^) 



It is possible to completely eliminate the r-dependence from this equation. (This is exactly 
as it should be since in the exact formulation, with all degrees of freedom retained, ^^i, 
and r have no counterparts.) To see this, we introduce the metric on the unit 4-sphere, 
'gjll^ = ^^j, (5*^(1)) , so that we may write 

?,.(S^(f))=f2?«. (3.15) 

Note that ^jil) is dimensionless, while ^^jy has the mass dimension —2. Writing D*^^-* and 
'g^^^ for the Laplacian and the determinant corresponding to 'g^u' we have, respectively. 



□ = r ^ nCi)^ and a/^ = r''^ ■ When we refer all quantities to the unit metric the 



average action reads 



47r G, 



+ FkirxB + rf[ 



(3.16) 



We may use this form of the functional on the LHS of eq. (13.141) . In the same notation, 
its RHS assumes the form 



Tr 



_□(!) 



X 



rxB)^fc2i?W(-(^^)+i^,'(fXB + f7: 
G,kdJG^HrxBreR^'\-j^) 



(3.17) 
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Notice that Xb, f, and r, both in (13. 160 and (13.171) . appear only in the cx-invariant com- 
binations Fxb and r /. It is therefore natural to introduce new fields, in terms of which 
r disappears completely: 

XB = rxB, l = rj, (f=r(j). (3.18) 

In analogy with </> = Xb + / we defined (f) = Xb + f- The new variables have the dimension 
of a length, while the old ones were dimensionless. The interpretation of xb, say, is clear: 
From g^j^ = Xb 9fj.u it follows that 

V = (fXB? gl^J = XB gli'J. (3.19) 

Hence xb is nothing but the radius of the sphere with the metric 'g^^. Likewise is the 
radius of a 4-sphere with metric ^^jy. 

To summarize: We have seen that we can eliminate the reference metric from the fiow 
equation determining the RG fiow of Gk and Fk. Therefore the corresponding beta func- 
tions cannot depend on it. Since the new variables are singlets under a-transformations, 
the fiow equation is invariant under the split transformations which express the arbitrari- 
ness of the reference metric. 

Having convinced ourselves that the fiow equations have the correct invariance prop- 
erties we shall now simplify our notation and omit the tilde from the new fields or, what is 
the same, we shall stick to the old fields but set r = 1 in all formulas. Then 'g^j^^ = "gjiJ and 
□ = D^i) so that we can omit the superscript "(1)" henceforth. The table of canonical 
dimensions changes accordingly; we have, along with [x'^] = 0, 

[l^iv] = [g,,u] = [ v] = [%u] = 

[/] = [7] = [xb] = [0] = -1 (3.20) 
[Hk] = [Fk] = -2. 
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With these simphfications the flow equation which we have to evaluate further reads 

(-□ + xl k' R''' (-^) + (xb + 7 (x)) ) " (3.21) 

x{(1-|^n) i?(0)(--Q,)-(--g^)i?(0)'(_^)- 

Here we used the anomalous dimension related to Newton's constant, defined as in [7], 

r]^ = kdk\nGk. (3.22) 

We shall see that the beta functions depend on the topology of the background 
manifold. We shall analyze the cases and in turn. For S^, we derived eq. (13.211) 
where ^^i, refers to the unit sphere now. For it is easy to see that the corresponding 
flow equations also can be obtained from (I3.2ip but with interpreted as the flat metric 

3.3 The RG Equation for Fk 

3.3.1 The 5^ Topology 

We begin by deriving the flow equation for the potential Fk{(f)). As explained above, 
we evaluate (13.211) for / = to this end, retain the full dependence on the (constant) fleld 
Xb = 4>, however. (Recall that = xb + / so that = Xb now.) We perform the trace 
on the RHS of (13.211) in the eigenbasis of □, the Laplace-Beltrami operator on the unit 
S'^. Its eigenvalues —£n and their degeneracies Dn, n = 0, 1, 2, ■ ■ ■ , are given by [60] 

Sn = n{n + 3), Dn = Un + l){n + 2){2n + 3). (3.23) 



With cr^ = J d^xy^ = 8n'^/3 the volume of >S"^(1), this leads us to 

3 a 4 



(A;9,F,(0) - r^N Ffc(0)) = Ti(0) + T2(0) (3.24) 
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with the two spectral sums 



T2 = --7]t^ip 2^ A 



(3.25) 



Here we employed the abbreviation 

ip = k(f). 



(3.26) 



Note that (f is dimensionless as has the dimension of a length. 

In principle fl3.24p gives the beta function of Fk for an arbitrary shape function 
R^'^\ In order to simplify the eigenvalue sums it is convenient to use the optimized shape 
function [61], however, the same that was also used in [I]: 



(l-z) 9{l~z). 



(3.27) 



With this choice, the final result for the flow equation of the potential reads 
^ f., ^ f., Gk (1-W2) P0V(A:0) + lr/Np(fc0) 

kd, F,(0) = F,(0) - wwTTm • 

Here p and p are two new, simpler spectral functions, for an arbitrary ^^i, defined by 



(3.28) 



p(V9) = Tr [0(^2 ^ Q^J ^ = Tr [(-□) 9{^' + □) 

In the case at hand, for the sphere, 

00 00 



n=0 



n=0 



These functions are of the form 



p(V5) = ^4('^max(V5)), 



(3.29) 



(3.30) 



(3.31) 



By definition, nmax(v^) is the largest integer n = 0, 1, 2, ■ ■ ■ such that £n = n{n + 3) < y?^. 
Here we introduced the finite sums 

N N 

J^{N) = J2Dn, MN) = J2^nDn. (3.32) 



n=0 



n=0 



18 



They can be worked out explicitly: 

j^(N) = ^N^ + lN^ + §N^ + lN + l (3.33a) 

J^{N) = -L + I + fl AT^ + f + f + I AT. (3.33b) 

Since the largest contributing eigenvalue has the quantum number 

{Lf for (/? ^ 1 
(3.34) 
for v9 < 1 

the results (13.33^ imply the limits 

^09^^ for 09 > 1 

p(^) ^{''^ ^ (3.35a) 

1 for (y9 -C 1 

,j^(p^ for > 1 
p{^) ^{''"^ ^ (3.35b) 

for <^ 1 
In particular we see that 

p(0) = 1, p(0) = 0. (3.36) 

On 5*^ the spectrum of □ is completely discrete. As a consequence, p{(p) and p{(p) 
have discontinuities at the y9-values at which nmax(v^) jumps. When one analyzes the 
differential equation 03.281) smooth functions p{ip) and p{(p) clearly would be easier to 
deal with. For this reason we henceforth adopt a certain smoothing procedure which 
replaces the original functions p and p by smooth interpolating functions. We shall again 
denote them by p and p. The way how this interpolation or smoothing is done is in no 
way unique. Choosing a specific smoothing procedure has the same conceptual status 
as choosing a particular cutoff function R^^^: It amounts to specifying how precisely the 
transition from the high to the low momentum regime takes place, i. e. how the modes 
start getting suppressed when their eigenvalue passes the threshold value given by k. 
Observable quantities derived from the RG flow must be independent of both R^^^ and 
the smoothing procedure. 
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It is convenient to assume the smooth p and p of polynomial form. We get a function 
with the correct behavior for ip and ^ 1 if 

4 

p((/?) = Ofc (/j'^, with ao = 1, a4 = 1/12 (3.37) 

fc=0 

and similarly for p. A good fit to the p{ip) of fl3.3ip is obtained with ai ~ —0.0118142, 
02 ~ —0.0832909, 03 ~ —0.000333389. For many purposes the following approximation 
is perfectly sufficient: 

p(y,) = l + 3^^^ p(^) = i^6_ (3 38) 

In particular when one searches for fixed points the fiow equation is needed in 
dimensionless form. For this reason we represent as 

Fk{<P) = k-^ Yk{k<P). (3.39) 
Here is a dimensionless function of a dimensionless argument /c = v?, i. e. 

Yk{^) = k''Fk{^/k). (3.40) 
The RG equation for Yk follows from (13.281) : 

kdkYk{(p) = {2 + 7]^)Yk{ip)-ifYi^{ip)~ — ^2 ^ Y^'iif) ■ ^ ' 

Here Qk = k'^ Gk is the dimensionless Newton constant. In terms of Yk, the original 
potential is given by Uk{4>) = Yk{k(f)). 

3.3.2 The Topology 

If the manifold on which the metrics ^^i, and 'g^^ are defined has the topology of 
a plane we may use eq. (I3.2ip as well, but with ^^i, = The operator □ is easily 
diagonalized in a plane wave basis then, and the spectral functions (I3.29P are found to be 

P(^)= (^) (/ d'x)/a, (3.42a) 
pM= (^) (/ d'x)/a,. (3.42b) 
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The resulting flow equation for the dimensionful potential reads 



Gk 



(I-I^n) 



06 



(3.43) 



247r 



02 ^2 + • 



For the dimensionless one finds the completely explicit partial differential equation 



kduY,{^) = {2 + r^^)Yu{v)-^Yl{^)-^ {l - \m) o (3-44) 



The equations (13.431) and (13.441) coincide exactly with the corresponding ones 
for 5*^, eqs. (I3.28P and ( I3.4ip . if in the latter the ^ 1- approximations from (]3.35ap and 



(I3.35bp are useqjfor all values of ip. This was to be expected, of course, because for ^ 1 
the spectral sums of S'^ are dominated by very many densely spaced eigenvalues which 
form a quasi-continuum. The difference between i?^ and S"^ is most pronounced for small 
values of (p. Here the finite volume of and the resulting discreteness of the spectrum 
strongly modifies the spectral density p. In the continuum it is proportional to v?^, while 
on the sphere it approaches p(0) = 1 for ip> <^1. Note that ip = k(j) can be interpreted as 
the radius of the sphere with metric measured in cutoff units. For ip <^ 1 the trace 
on the RHS of the FRGE is dominated by a single eigenvalue, namely the zero mode of 
□ which has n = 0. For the topology, on the other hand, the p ~ (/^^-behavior extends 
down to = 0. 

In [ I ] we explained in detail how the special status of the metric affects the RG equa- 
tions and why those for conformally reduced gravity are different from the corresponding 
flow equations in scalar matter fleld theories. It is therefore instructive to compare the 
partial differential equations (13.431) and (I3.44p to their scalar counterparts. In a standard 
scalar theory, the analog of eq. (13.431) for F^, say, would have the last factor on its RHS 
replaced by k^ / [k"^ + F^'((/))). It does not contain the crucial extra factors of (f) which 
originate from the xs-dependence of TZk- It is clear that these extra factors of modify 
the RG evolution of the function Fk{4>) quite signiflcantly. 



■*Note that the first factor on the RHS of (|3.42ap and (j3.42bp . respectively, coincides exactly with the 
S'^ results (|3.35a[) and (j3.35bp valid for ^ 1. The second factor / d^x/a^ takes the different volume 
normalizations in the two cases into account. 
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Within the present framework, the CREH approximation of [I] amounts to the 
following more restrictive ansatz for Y^: 



Here 



for i?^ 



Co 



(3.45) 



(3.46) 



1 for 

and Afc = Ak/k"^ is the dimensionless cosmological constant. We can derive its RG equation 

A;9,Afc = /3,^^™(^7fc,Afc) (3.47) 

by inserting fl3.45p into fl3.4ip or fl3.44p and comparing the coefficients of ip^. From either 
equation one obtains 

a,) = -(2 - r^^) A, + £ (1 - I r^N) ^3^. 
This is exactly the beta function found in [I]. 



(3.48) 



3.4 The RG Equation for Gk 

In order to close the system of equations we need an equation for dk Qk or, equiva- 
lently, the anomalous dimension: 



kdkgk= 2 + ri^{gk,[Yk]) 



9k- 



(3.49) 



In the case at hand, r/N is a function of Qk and a functional of Y^- The desired expression 
for 1]-^ is obtained from eq. fl3.2ip by fixing a constant value of xb and performing a 
derivative expansion of the functional trace up to order (df)^. 



3.4.1 The Topology 

For ^fj^i, = 6f^u the result for r/N can be inferred from the calculation in Appendix A 
of [I]. Expressed in dimensionful quantities, it reads 



Gk 

ovr 



k(f)i 



P 02 



1 + 



Gk 
12^ 



kcpi 



(3.50) 
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Here S4 and S4 are threshold functions defined in [I]. They depend on R^'^\ For the 
optimized shape function fl3.27p one has, for instance, 



1 



1 



0. 



(3.51) 



4 (1 + 

In (I3.50p the derivatives of Fk are evaluated at a fixed field value, 0i. It is the value of 
Xb about which we perform the expansion of fl3.2ip . There is a certain arbitrariness as 
to the choice of xb = 0i in '7n- It is due to the fact that the truncation used allows 
for a field independent wave function renormalization Z^k oc 1/Gk only, see Appendix A 
of [I]. In standard scalar calculations (pi is usually identified with the minimum of the 
potential [33]. Notice that since 0i is fixed, ?7n is indeed a functional of and not a 
function of the field. 

With the optimized shape function rj-^ assumes the simple form 

2 
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1 + 



F!'( 



P 02 



(3.52) 



We see that r/N is negative if Gk > 0, for any choice of (pi and for any potential function 
Fk{ ■ ). This means that the gravitational interaction is aniz screening: As we increase k, 
Newton's constant decreases according to kdk G^ = r^N Gk- 

In terms of the dimensionless quantities Qk and the anomalous dimension (13.521) 
reads, with ipi = k(pi, 
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If we insert the CREH form of the running potential, V^creh 



from (13331) : 



2 

'3^ (T 



2 A. 



(3.53) 

{if) = — |AfcV5^, we obtain 
(3.54) 



This is precisely the anomalous dimension of the CREH approximation which was derived 
in [I]. (In [I] the above r/^'^™ was denoted Note that (13.541) is independent of the 

expansion point (fi. exactly when Yk oc it drops out from the formula for t^n- 
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When we solve the coupled flow equations later on we set y^i — oo for the expansion 
point. The motivation is that the relevant solutions Yk{(p) turn out to behave as oc ip^ 
for if —>■ oo. Hence with this choice % is given by the CREH formula fl3.54p which is 
independent of the precise value of (fi as long as it is sufficiently large. Other choices of 
(pi are conceivable, for instance minima of the potential. It turns out, however, that at 
least in the particular interesting UV fixed point regime Yk has no stationary points at 
all. 

3.4.2 The 5^ Topology 

Also in the S'^ case we shall send ipi to infinity. The corresponding sphere has a 
very large radius then and we can expect that t^n is well approximated by the result. 
For this reason we are going to employ fl3.54p also in the spherical case. We defer a more 
precise treatment, including an evaluation of t^n on a sphere of finite radius, to a future 
publication. 



4 General Properties of the LPA Flow 

The RG fiow on the infinite dimensional theory space with coordinates [g,Y{-)) is 
described by the coupled equations 

kdkYkiip) = PY{9k,Ykiip)) (4.1a) 
kdkgk = Pg{gkA^k]) = 2 + r]t^{gk,[Yk]) gu- (4.1b) 

The beta function for the potential, /3y, is given by the RHS of fl3.4ip for S'^ and fl3.44p for 
R^. The anomalous dimension will be used in the form fl3.53p . Its functional dependence 
on the potential is in terms of derivatives of Yk evaluated at the fixed expansion point ipi 
which we shall leave unspecified for the time being. 

Alternatively we may use the fiow equations in dimensionful form. The partial 
differential equation for -Ffc(0) = k^'^Yk{k(f)) is given in fl3.28p for S^, and in (13.430 for 
R^; in either case Gk evolves according to kdk Gk = Tj^Gk with the r^N of (13.521) . 
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In this section we are going to discuss various general properties of the coupled RG 
equations. 



4.1 The — ^ oo Asymptotics 



Let us separate off the CREH form of the potential, 
from the function Y^: 

Y,icp)^AY,i^) + Yr-i^). 



CREH 



{(f) = Co (y?^ - 1 Afc ip^ 



(4.2) 



Exploiting (I3.47P we can derive a flow equation for the "correction term" AYJj. For we 
obtain 



kdk AYkif) - (2 + t/n) AYk{(f) + (f AY^{<f) 
Here we abbreviated 



(4.3a) 



Bfc(<^) = 



12 



Vn 



_pM\ 



l-2Afc+ 



2 + An"(v^) 



2 36 



(4.3b) 



(1-2 A, 



What prevents us from setting AY^ = is the inhomogeneity on the RHS of eq. fl4.3ap . It 
controls the type of terms which are generated during the RG running. It is important to 
notice that, in the limit ^ oo, this inhomogeneity simplifies under certain conditions. 
In fact, from fl3.35p we know that p V^'^/12 and (p — 'p/^p^) — > (p^/36 when (p oo. 
Therefore we learn from (I4.3bp that Bfc((y9) vanishes for ^ oo if AYl!{ip) / ip^ for 
if oo. The latter condition means that the asymptotic growth of AYi^{ip) is weaker than 
oc v?^. Let us assume that at some initial value of k the potential AYfc((y9) does indeed 
grow more slowly than ip^. Then the RHS of fl4.3ap . for large values of (p, is proportional 
to (pP' and therefore does not generate any components (monomials) in AY^ that would 
grow oc or faster. This implies that if initially A Y^ {p>) contains no terms growing like 
p)^ or faster then no such terms will be generated by the flow. 

While eqs. (14. 3 p were written down for it is obvious that this statement holds for 
both and R\ 
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If AYk grows more slowly than (p^ then, at any k, the large-y? asymptotics of the full 
potential is always given by the CREH potential or, more specifically, by its cosmological 
constant term since the v3^-piece is subdominant: 

Mv) -| Afc for ip^(X). (4.4) 

4.2 Structure of the Initial Value Problem 

In the case at hand the specification of the truncated theory space involves choosing 
an appropriate space of functions on which evolves, in particular boundary conditions 
must be defined for Y^. The above discussion shows that it is consistent with the differ- 
ential equation to partially define the theory space by the requirement that Yk{(p) does 
not grow faster than (f'^ for y9 — oo. In the following we shall adopt this choice. It is 
further motivated by the RG fixed points we shall find in the next section. In particular 
a NGFP suitable for the asymptotic safety construction and all corresponding "physical" 
trajectories (those hitting it for k oo) are inside the space of functions growing not 
stronger than (p^. 

As we mentioned already, the asymptotics (14.41) motivates the choice (fi-^oo for 
the expansion point in the calculation of r/N- We shall adopt this choice in the following. 
It has the technical advantage that now r^N is no longer the complicated functional of Y^ 
given in (13.531) . but rather assumes the simple CREH form (I3.54p . It depends on gi^ and 
only one particular characteristic of Y^, namely the coefficient appearing in its large-*/? 
asymptotics (14. 4p . 

With the above two choices made, solving the initial value problem for the system 
(I4.1b .b) proceeds as follows: 

(i) Fix, at the initial scale fcinti, a Newton constant (yf^j^^y and a potential Y^.^^^^^ip). 

(ii) Extract the number A^.^^y = —6 lim<^^oo ^fcinti if) I f'^ from the asymptotic behavior 
of the initial potential. 
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(iii) Solve the two coupled ordinary differential equations of the CREH truncation 

kd, = pr^"" (^^r™, a^^™) 



CREH 
Ok 



(4.5) 



subject to the initial conditions g^^^^ = Qh^ui ^kf™ = ^huu- ■'■^ (14.51) . r]^™ and 
^CREH g^j^g given by eqs. (I3.54p and (13 ■481) . respectively. 

(iv) Use the solution of (14. 5 p to express the anomalous dimension as the following explicit 
function of k: 

^N(A;) = r^r"(^r^Ar"). (4.6) 

(v) Insert the solution = g^^^^ of (14.50 and the anomalous dimension ?7n = ri-t<{{k) of 
]6]) into the partial differential equation for Y^, eq. (I3.4ip . and solve it with the 



initial condition specified. 

Note that by this procedure the equations for g^ and get decoupled; at step (v) we are 
dealing with an equation for alone, albeit one which contains complicated /c-dependent 
coefficient functions. 

4.3 A Robustness Property of the CREH Approximation 

The structure of the initial value problem for the LPA of conformally reduced gravity 
sheds an interesting light on the simpler CREH truncation. Both of them retain only the 
simplest derivative term, {dipY- The former allows for an arbitrary potential Yk{Lp), the 
latter only for the y9^-monomial. In this sense the LPA represents a refinement of the 
CREH approximation. 

Usually when one refines a truncation ansatz by adding further field monomials to 
it the projection of some RG trajectory "living" in the new, higher dimensional space 
onto the smaller space of the old truncation equals only approximately (or in the case of 
badly chosen spaces, not at all) an old trajectory in the smaller space, as computed with 
the simpler beta functions of the lower dimensional flow. Stated the other way around. 
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if we compute the flow of a smaU set of couplings in a simple truncation, then reflne the 
truncation, and flnally project the higher dimensional trajectories thus obtained onto the 
old, smaller space, we flnd that the projections of the higher dimensional trajectories will 
in general not coincide with the lower dimensional ones. (But they will approximately if 
the truncations are reliable.) 

The above discussion of the initial value problem shows that the LPA of the con- 
formal factor, seen as a reflnement of the CREH approximation, is a notable exception 
to this rule: The projection of the RG trajectories in the inflnite dimensional (^g,Y{-))- 
space onto the 2-dimensional {g, A)-space coincides exactly with the trajectories computed 
from the two beta functions (3g and Px appropriate for the smaller space. Generalizing 
the space of potential functions from the 1-dimensional line {A<^^} to the inflnite dimen- 
sional {Y{ ■ )} has no impact on the RG evolution of Qk and A^. This indicates once more 
that the Einstein-Hilbert truncation is particularly robust under the inclusion of further 
invariants. 

This robustness property has, among others, the following implication. We know [10] 
that some trajectories of the Einstein-Hilbert truncation, those of Type Ilia, terminate at 
a flnite scale, and the same was found to be true for the corresponding CREH trajectories 
[I]. Because the LPA trajectories [gk,Yk{-)) are generalizations of CREH trajectories 
idk, Afc) in the sense explained above, the breakdown of a CREH trajectory at some fcterm 
implies that the corresponding LPA trajectory terminates at the same 
Generalizing the (f^-teim to an arbitrary function Yk{ip) seems not to help in extending 
the validity of the truncation towards the IR. This result is somewhat surprising since 
in the case of scalar matter fleld theories trajectories which break down at a flnite scale 
within a polynomial truncation (those in the broken phase) typically can reach k — 
when the LPA is employed. 

4.4 The cp ^ Behavior 

Next we analyze the partial differential equation for Yk{(f) in the limit </? — > 0. We 
flrst focus on the inhomogeneous term of the flow equation, the last term on the RHS of 
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fl3.4ip . We may assume that for ip ^ the function p{<f) vanishes at least proportional to 
Hence, in the numerator of the inhomogeneous term, (1 — ip^ p{ip) + | r/N p{v) = 

(9(<y5^). In order to estimate the denominator ip"^ + Y^'{if) we decompose Yk according 
to Yk{ip) = Y^^^{{p) + Y'^^"^((y9) where the "regular" part Yj^'^^ is defined by the condition 
{Y^^^)"{(f) < oo for y9 — s> 0, while for the "singular" part V^"^ the second derivative blows 
up for vanishing field: (Y^^™^) oo for (f ^ 0. Let us assume a behavior of the form 
{Y^^'^y'iip) = a/ip^' with constants a and /i > 0. Since, by definition, (F^^^)" is 
or vanishes even this implies the following structure of the inhomogeneous term: 

?M = = oi^'^n- 

We conclude that the inhomogeneity in the flow equation vanishes oc v^^'*'^, i.e. with 
a power larger than 2, when 0. Therefore, in a small-y? expansion it contains only 
terms whose second (/^-derivative vanishes for — > at least like y?^. Hence, upon inserting 
Yk = Y^'^^ + y^^"^ into (13.411) . we see that there are no terms in the inhomogeneity that 
would match those in Yj^^^^. As a result, the "sing" part of the potential has no source 
term; its RG equation is homogeneous, implying that no "sing" terms are generated if 
they are absent originally: 

¥^^" = =^ kdkYr^ = Q. (4.7) 

Note that the attributes "regular" and "singular" refer to the second (/^-derivative of 
Yki not the potential itself. As an example consider Y{ip) = ip"^ with an exponent m G R. 
Then Y"{(p) = m (m — 1) (p^'"^ so that Y is "regular" if m = 0, or, m = 1, or m > 2; 
otherwise it is "singular". 

An argument similar to the above one has also been used in [21], in a different 
physical regime though (i? 0, corresponding to ip ^ oo here). 

5 The Fixed Points 

In this section we continue the analysis of the RG equations (14.11) by searching for 
fixed points (^g^,Y^{-)), i.e. solutions of /5y (^f*, (</))) = = (3g[g^, [Y*])- Depending on 



29 



whether Pg = is achieved by (7* = or r]t^{g*, [Y.^]) = —2 we call them, in a slight abuse 
of language, a Gaussian fixed point or a non-Gaussian fixed point, respectively. (For a 
more precise discussion of the distinction see [29].) 

5.1 The Gaussian Fixed Point (i?^ and 5^) 

We look for fixed points with = 0. As a consequence, the second condition /3y = 
boils down to 

(2 + r/,)n(^)-^i;'(^) = 0. (5.1) 

This equation obtains from (13.411) as well as from (13.441) . so it holds for both and Rf^. 
Here i]^ = ^'^{g*, [Y*]) is the anomalous dimension at the fixed point. By eq. (13.531) it 
vanishes: rj^, = 0. The differential equation (15. ip is trivial to solve then, with the result 

Yn^) = c^\ gr = o. (5.2) 

Here c is an arbitrary constant of integration. 

As c in not fixed by the equations it might seem that we found a 1-parameter family 
of fixed points. However, let us recall the CREH form of the potential: Y^^^^{(f) = 
Co — (Afc/G) ip'^. Its quadratic term cq (p^ is linked to the kinetic term by a local Weyl 
rescaling, a "a-transformation" . (The cq (f'^-teim follows from the | R-teim in (13. ip upon 
inserting ^^i,; both the | i?- and the D-term emerge from J y/g' R when the conformal 
factor is separated off.) If we impose the condition of local a-invariance we must pick the 
value c = Co therefore; then can be identified with (a part of) the conformal reduction 
of J R- Thus c = for R^ and c = 1 for 5*^. It should be noted here that local 
(T-invariance really amounts to an additional assumption. After insisting on a ^^j^-metric 
of the S^-type the split symmetry was broken down to global cr-transformations, and they 
are not strong enough to forbid the generation of (/j^-terms unrelated to the kinetic term 
and / y/g' R. For conceptual clarity we shall leave c arbitrary in the following, but clearly 
c = Cq is the most sensible choice. 
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5.2 Linear Stability Analysis 

Given an arbitrary fixed point )), not necessarily a Gaussian one, we can 

explore its stability properties by linearizing the flow in its vicinity according to 

gk = g* + Sgk (5.3a) 

Yk{^) = Y,{^)+6Yk{f). (5.3b) 
For the "tangent vector" [6gk, SYk{ ■ )) we make the ansatz 

Sgk = £y9i-^j =eygQ~^* (5.4a) 

5Y,{^) = e T{^) (jj T(y,) e-^* (5.4b) 

Here ko is an arbitrary fixed scale, t = ln(A;//co) denotes the "RG time", and e is an 
infinitesimal parameter. We shall determine [Sgk, SYi^^ ■ )) by expanding the flow equation 
to first order in e. Then (y^, T( ■ )) is an eigenvector ("scaling field") of the corresponding 
infinite dimensional stability matrix with the negative eigenvalue ("critical exponent") 6. 
Scaling fields with Re^ > (Re^ < 0) are said to be relevant (irrelevant). They grow 
(are damped) along the RG flow, i. e. when k is lowered. 

By definition, the fixed point's UV critical manifold S^uy, or synonymously, its 
unstable manifold, consists of all points ((?, Y{ ■ )) which are pulled into the fixed point by 
the inverse fiow. The tangent space to ^uv at the fixed point is spanned by the relevant 
scahng fields. Hence Auv = dim^uv equals the number of relevant scaling fields, i.e. 
the number of positive Re^'s. (See [29] for a more precise discussion.) 

The scaling fields and dimensions are to be determined from the following two equa- 
tions which are obtained by varying (14. lap and (]4.1bp . respectively. The first one is 

[2 + 9 + T]^] 6gk + g* 5?7n = (5.5a) 
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and the second equation reads 

^ 5gk (1 -^*/2) v?V(<^) + lv*p{f) 



'47r '''' v' + Y:'{^) 

27r b2 + y;%)f 



(5.5b) 



Here 77* = ?7n(5'*, [K*]) and ^r^N = m{g* + ^fiffc, + ^Ykl) - V*- As is stands, the last 
formula holds for S^. The one for R"^ follows by the replacement 

pif) — ^^</'^ pif) — 'Tsf^- (5-6) 

Next we apply this analysis to the GFP. 

5.3 Stability Analysis of the GFP (R^ and 5^) 

For the anomalous dimension given in (I3.53P we obtain at the GFP fl5.2p : 

^N(^r, [rr]) = 0, {Srm){9r, [Yrd = o. (5.7) 

This turns f l5.5ap and fl5.5bp into the following two conditions for (y^, T( ■ )) and 9: 



[2 + 9]y, = (5.8a) 

[2 + e]T{^)-^n^) = I^^^P^ (5.8b) 

The solutions to this linear system are easily found. As f l5.8ap can be satisfied by either 
T/g = or 2+0 = we have two types of scaling fields at the GFP. They can be summarized 
as follows: 

(i) For every 9 E H there exists a scaling field {yg = 0, T = Tg) with scaling dimension 
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9 ancj] 



(5.9) 



with constants cq. For 9 > 0, 9 = 0, and 9 < these scahng fields are relevant, 
marginal, and irrelevant, respectively. 

(ii) There exists a single additional scaling field {yg 7^ 0, T = T_2) with scaling dimen- 
sion 9 = —2 and 



This scaling field is irrelevant. 

The scaling fields of type (i) are remarkable in that a monomial n G R, has the 
dimension 



Relative to a standard scalar field theory, the scaling dimension at the GFP is shifted by 2 

units. This shift can be traced back to the additional factors of (p in the FRGE which, in 

turn, originate from the xs-dependence of TZk which was needed in order to give the desired 

physical interpretation to the cutoff and to implement "background independence" . Again 

we see that the KG behavior of conformally reduced gravity, despite its appearance, is 

very different from that of a scalar matter field theory. 

In Table [U we list the scaling fields with integer exponents n as an example, along 

with their dimensions. The two entries printed bold-faced are the scaling fields accessible 

to the conformally reduced Einstein-Hilbert truncation. The scaling field {yg = 0, (p^) 

with 9 = +2 corresponds to the 2-component eigenvector of the stability matrix for the 

{g, A)-plane which is parallel to the A-axis; the other one, {yg 7^ 0, T_2) with 9 = —2, 

represents the eigenvector with both a non- vanishing g- and A-component. This pattern 

^Here we assume that the exponent n oiTe (x is reaL At this point it is not clear though whether 
one should impose 7i e R or Z or M, or one should allow all n G C even. We are going to discuss this 
issue of the proper choice for the {Y{ ■ )}-function space in detail in Subsection 15.51 




(5.10) 



9 = n-2. 



(5.11) 



33 



e 




-3 


-2 


-1 





+ 1 


+2 


+3 






irrelevant 


marg. 


relevant 


yg = o, T = 
























T_2 





Table 1: Scaling fields at the Gaussian fixed point 

is exactly the same as in the full [9,10] and the conformally reduced [I] Einstein-Hilbert 
truncation. In fact, T_2 is essentially oc (/?^, plus corrections which are subleading at large 
if. If we perform the integral (15.101) with the spectral function p of (15. 6p . appropriate for 
the case, we obtain 



T-2(^) 



967r L 



(5.12) 



If we invoke local cr-invariance and set c = Co[-R'^] = we find that T_2 oc ip^ on i?^. 

On the smooth p{(f) given by the polynomial (13.371) yields the analogous scaling 
field (if c > 0): 



Ztc 



1 
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TE'-P + —'^ + 2 1*^2-^) +(ai-2ca3)v9 



(5.13) 



+v2c (2 caa — ai) arctan(y9/v2c ) 



+ ^ ri-2ca2 + y) ln(v9^ + 2 c) + const 



We see that T_2 becomes a pure yj'^-term for ip ^ 1. For </) — > there would seem to 
be a singularity in the special case c = 0. But fortunately the distinguished value is 
c = co[S'^] = 1, and for this value T_2 is regular for all > 0. For it approaches a 
constant. 

5.4 The non-Gaussian Fixed Point 

Let us now try to find fixed points {g^,Y^) with rj^, = —2. We begin with the flat 



case. 
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5.4.1 The Topology 



According to fl3.44p and fl3.53p the coupled system to be solved is, with ipi fixed but 
unspecified for the time being, 

-,6 



9* 



487r 



1. 



(5.14a) 



(5.14b) 



In order to analyze fl5.14al) it is convenient to introduce h{Lp) = YI{lp) which satisfies 
the first order equation 

n5 



9* 



(5.15) 



.ISvr/ v?2 + /i'(v?) 

An asymptotic analysis of this equation reveals that, to leading order, its solution has 
essentially the same behavior for (/? ^ 1 and ^ 1: 



-\L^^ if(/?>l 



(5.16) 



Here L and L are constants which must satisfy the same quadratic equation. 



l2nL (1 -2L) = g,, 



(5.17) 



and likewise for L. This equation can have two different real solutions, so L and L can 
be different in principle. The behavior fl5.16p motivates introducing a function W{ip) for 
< (/? < oo by 

2 



The new function satisfies the boundary conditions 



iy(o) = 1, 



W{oo) = L/L. 



(5.19) 



35 



Inserting flS.lSp into flS.lSp we find that the relevant differential equation has the structure 

(p-^W{ip) = B{W{ip)) (5.20) 

with the "beta function" 

B{W) = 2 LW - {1-2 L) /w]. (5.21) 

2 L \- J 

Being first order in ip, the equation (15.201) with both boundary conditions (I5.18P imposed 

is overdetermined and the existence of a solution is questionable a priori. If we start the 

integration of (I5.20p at = with the initial condition W{0) = 1 then the resulting 

solution W{(p) has no reason to approach L/L for —>■ oo, at least for a generic B{W). 

However, the function (I5.2ip has a special property which implies the existence of a 

solution. It has zeros at W = 1 and W = L/L: 

B{1) = 0, B{L/L) = 0. 

(To verify the second zero one must exploit that L and L satisfy (15.171) .) Obviously 
= 1 is a "fixed point" of the "flow equation" (I5.20p : If we impose the initial condition 
1^(0) = 1 and integrate towards larger ip^s we obviously get 

W{ip) = 1 for < < oo. (5.22) 

This solution also satisfies the second boundary condition W{oo) = L/L ii L and L are 
the same solution of the above quadratic equation. (Starting from the other "fixed point" 
W = L/L and integrating backward from ip = oo leads to the same conclusion.) 

Thus, integrating (I5.18P with (15.220 . we have shown that the fixed point potential 
is of the form 

Y,i^)=y,-lL^'' (5.23) 
where L is a solution of (I5.17p . and ?/* is a constant. 



The second condition {g*,Y^,) must satisfy is fl5.14bp . Inserting the result (I5.23P for 
Y*((/)) it boils down to 
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Note that, thanks to the specific form of y*, the expansion point ipi dropped out of this 
equation. 

What remains to be done is to solve the two coupled algebraic equations (15.171) and 
fl5.24p for the constants and L. Actually, with the identification L = A*, these two 
equations are exactly the same as those we encountered in Subsection 5.2 of [I] where 
we determined the NGFP (5'*,A=k) within the CREH approximation (using the "kinetic" 
version of the anomalous dimension, ■'7^™''). As a result, the much more general LPA 
yields the same fixed point values for g and A as the conformally reduced Einstein- 
Hilbert truncation: g^, = g^^™, L = A* = A^^"^". Recalling the values found in [I] we can 
summarize our result for the NGFP as follows: 

y*{v) = y* - 7: K '^'^ (5.25a) 
o 

A. = I j^^^^ ^ 0.279 (5.25b) 

91/3 

= 67r ^ ^ 4.650 (5.25c) 

(1 + 21/3)2 

The constant is not determined by the flow equation. Except for this constant, the 
equations have not taken advantage of the possibility to generalize the fixed point potential 
beyond the functional form of the CREH approximation. This result is quite remarkable. 
It might be related to the impressive robustness and stability properties the full Einstein- 
Hilbert truncation is known to possess [9, 11]. 

5.4.2 The 5^ Topology 



In the case of the condition ?7n = — 2, again, translates to eq. (15.14bp . while the 
other fixed point condition /5y = 0, by fl3.4ip . becomes 



with the effective p-function 



l':(^) = -^4f#T (5.26) 



Peff(¥^) = P(V^) - (5.27) 

2 
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For — > oo, pcfr approaches </''*/18 so that (15.260 coincides with its counterpart, eq. 
f l5.14al) . in this hmit. This imphes that the asymptotic form of the solution to (15.261) 
agrees with that of the solution: 

y.(y^) > -lK^' = Y^yrni^). (5.28) 

Inserting = Fasym into (15.261) and letting ^ oo we obtain a first relation among the 
constants and A,,; it is given by eq. (15.171) with L = 

Furthermore, the fixed point condition r/N = —2 has the explicit form ( I5.14bl) . We 
adopt the choice (fi-^oo here. As a result, we may replace by yasym in ( I5.14b[) . For 
(fi^oo the point (fi actually drops out as before, and what remains is a second relation 
among g^, and A*; it is given by eq. (I5.24p with L = A*. The conditions (15.171) and (15.241) 
determine g^, and A* uniquely. Therefore we can conclude that, exactly as for R^, the 
constants A* and g^: assume precisely their CREH values (I5.25bl) and (I5.25cl) . respectively. 

So, what is left to be done is to solve the ordinary differential equation (I5.26P for 
the by now known value of g* subject to initial conditions which involve A^=. We require 
that for large ip the function Y^ agrees with Fasym, and likewise their derivatives. We start 
the integration at some (p where we impose the initial conditions 

Yi^) = Fasyml^) ^'l^) = Ksym(^)- (5-29) 

From if we integrate backward towards smaller (p. Ultimately we are interested in the 
limit ^ ^ oo. 

The existence of a fixed point potential is by no means guaranteed. It could happen 
that the function one obtains by integrating (I5.26p . (I5.29P towards smaller ip develops an 
unacceptable singularity at some point. Then we would have to conclude that there is no 
fixed point for Y^. The potentially dangerous feature of the differential equation (I5.26P is 
the denominator (f"^ + Yl'{if) which could possible vanish at some ip. 

Let us assume this does not happen and there exists a solution Y^, which is well 
defined for < (p < Ip oo. Then this function is necessarily monotonically decreasing, 
i. e. Y^{ip) < everywhere. The reason is as follows. A non-singular solution has + 
Y"{ip) > everywhere, and also the function Pcs{<p) is found to be positive for any ip. 
Since g^, > 0, this entails that the RHS of the fixed point equation (15.260 is negative. 
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Unfortunately it is not possible to solve the differential equation for analytically 
in closed form. We shall have to resort to numerical techniques therefore. It is easy, 
however, to find the leading behavior for <^ 1. Making a power series ansatz and 
working out the coefficients one gets 

y*{^) = yo~yiv- + 0{^'). (5.30) 

The constants of integration yo and yi are undetermined at this stage; ultimately they 
must get fixed by matching (15.301) with the large-*/? solution. The monotonicity Y^{lp) < 
requires yi > 0. Eq. (I5.30p implies that Y"{{p) = 0{{p) which vanishes for — > 0. 

The numerical analysis of (I5.26P shows that there does indeed exist a solution which 
is regular for all > 0. It is displayed in Fig. [1] where it is compared to the CREH 
potential YJ~'^^^{(f) = ip^ — (A*/6)y9^ with the same value of A* (which is the same as 
^asym, of coursc). In the Figure we actually plot the negative fixed point potential — 
since, up to a positive factor, it agrees with the true potential [/*, the non-derivative term 
in die Euclidean F^. As expected, Y^.{ip) is indeed a monotonic function, contrary to the 
CREH potential. We also observe that is approximately linear for small y?, in accord 
with fICTD . 

It is important to observe that y*(v5) and Yl'{ip)) are finite for \ 0. The latter 
result means that the S'^ fixed point potential has a vanishing "sing" part in the sense of 
Subsection 14.41 The same is also true for as the explicit formula (15.23^ shows. 

5.5 Stability Analysis of the NGFP 

In this subsection we perform a linear stability analysis of the non-Gaussian fixed 
point. According to the asymptotic safety idea, the infinite cutoff limit of QEG is to 
be taken at this fixed point. As a consequence, the structure of the NGFP's ultraviolet 
critical manifold determines the physical properties of the theory in an essential way. 
In particular its dimensionality Auv = dim^uv equals the number of parameters it can 
depend on and which are not fixed by the requirement of a "safe" UV behavior. 
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Figure 1: The negative fixed point potential —Y^, for 5*^. The dashed hne is the CREH 
potential with the same value for A,.. The figure on the LHS is a detail of the one on the 
RHS; the latter shows that coincides with the CREH potential asymptotically. 



5.5.1 Scaling Fields and the Condition n' < 4 

We restrict ourselves to the topology here. Then the relevant fixed point solution 
is f l5.25p with the CREH values of and A*. The first condition for the scaling fields 
and dimensions is f l5.5ap . and the second one is ( ]5.5bl) with the spectral functions (15.61) 
inserted. 

Thanks to rj^: = —2 the first condition assumes the form 

{Sv^){9*, [y*\) = -0 6gk/g*. (5.31) 

The anomalous dimension is given in eq. ( 13.53^ . Varying it with respect to g and Y{-) 
while keeping (pi fixed we obtain after inserting the fixed point values: 

Using (I5.32P together with (15.40 in (I5.3ip we arrive at a condition relating the "compo- 
nents" of the scaling field, (y^, T( ■ )), and its dimension, 9. It is convenient to distinguish 



+ 



1 -2A. 



(5.32) 
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the cases i/g = and i/g 0. For i/g = the condition boils down to the ^-independent 
relation 



ipi 1 — 2 A* ifl 
For yg the analogous condition reads 



0. 



(5.33) 



9 = 2 



A* 



Svryg (1-2 A.) 



r"(v9i) , 8 A. T"(^i) 



1 - 2 A, 



(5.34) 



In the case at hand the second condition is found to be 



{-6) [T(^) - yg/g:\ + ^ T'(^) - a T"(^) = (71 ^ - 72) 



with the useful abbreviations 



a 



g* 



21/3 



187r (1 - 2 A*) 



0.41997- ■• 



71 = 



72 = 



A, 



1 



6 [g, 247r (1 -2 A, 
1 



(5.35) 

(5.36) 

(5.37) 
(5.38) 



187r (1-2 A,) 

In the special case yg = we get the homogeneous Eulerian differential equation 

a^^T"{^)-^T'{^) + eT{^) = 0. (5.39) 
We shall now analyze the cases yg = and ?/g 7^ in turn. 



(a) Scaling fields with tfg = 0. In order to solve the differential equation (15.391) we 
make a power law ansatz T{ip) oc y?", admitting arbitrary complex exponents n G C a 
priori. The ansatz is a solution if 6' = 6{n) where 



9{n) = {1 + a) n — an^. 
Solving for n we get the two possible values 

n±{e) =uj± Vcu^ - e/a 



(5.40) 



(5.41) 
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where 



l + a 
2a 



1.690 



(5.42) 



So, ii 6 ^ auP' , i. e. n_|_ ^ n„, the general solution of fl5.39p for a given d reads 



(5.43) 



where /i+ and /;._ are arbitrary complex constants. In the exceptional case Q = aiJ^ the 
general solution is instead {h\^i G C) 



(5.44) 



Generically both Q and n will be complex. We set Q = Q' ^ iO" and n = n' + in" 
with real and imaginary parts 9' , n' and 9" , n", respectively. Then f l5.40p decomposes as 



9' = (l + a) n' - an'^ + an"'^ 
9" = n" [l + a-2an']. 



(5.45a) 
(5.45b) 



We shall also need the inverse of these relations. After some algebra one obtains the fol- 
lowing results for and as functions of 9' and 9". Three cases are to be distinguished: 
(i) The case 6" 0: 



n'^i9',9")=uj 



1 ± 



V2 



9' 



+ 



+ 1 



l/2\ 



n'!,{9',9") = T^ sign(n 



9" 



9' 



1/2 



(5.46a) 



(iia) The case 0" = and 6' < a 



n'^{9\ 9") =uj 1 ± ^l~9'/auj^ 



nl[9\9") = Q 



(5.46b) 
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(iib) The case 0" = and 0' > a uj"^: 

n'^{e',e") = iu 

(5.46c) 

n'^{9', e") = ±uj ^/9'/auj^ - 1 

By now we imposed only one of the two conditions scahng fields must meet. The 
other one, for ?/p = 0, is eq. (I5.33p . If T{{p) oc it reads 

n (n - 1) [n - 2 + 12 a] ip^'^ = (5.47) 

where we used that 8A,,/(1 — 2 A*) = 12 a ^ 5.0395. Eq. fl5.47p is satisfied if at least 
one of the following conditions holds: n = 0, n = 1, n = 2 — 12 a, v?i~^ = 0. The first 
three of them provide us with finitely many scaling fields only, while we expect infinitely 
many, of course. Therefore we should demand that v'""^ = in an appropriate sense. 
This condition implies that either (fi = with Re n > 4, or that (fi^oo with Re n < 4. 
The first case is clearly excluded; as the metric with = is singular one should not use 
<y9i = as the expansion point in the computation of ?7n. It remains the option — oo 
which we had actually advocated for a different reason already: If ip is large enough, ?7n 
becomes actually independent of it. 

Thus we conclude on the basis of the stability analysis that our present approxima- 
tion can lead to a consistent picture only with the choice (fi-^oo. For the scaling fields 
with yg = 0, T{ip) oc this entails that the exponent is constrained by 

Ren{e) = n'{e',e") <4. (5.48) 

In order to find the spectrum of scaling dimensions we now combine the two condi- 
tions scaling fields must meet. The allowed values of {$', 9") are such that the real parts 
given by the expressions (15.461) are strictly smaller than 4. After some tedious algebra 
one finds that there exist the following two families of scaling fields (y^, T) = (0, </)"): 
The "n_|_-family": There exists a scaling field (?/g, T) = (0, (y9"++*"+) of complex scaling 
dimension 9 = 9' + i9" for every point (6*', 9") in the complex 6'-plane such that 

9' > -ai + a2 9"^ (5.49) 
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where ai =4 {2^/^ -I) ^ 1.0397 and 02 = 3 ■ 2^/3 (7- 2^/3 _ 3^-2 ^ 0.III6. The associated 
exponents n' + in" are given by the eqs. fl5.46p with the upper sign. The domain of allowed 
scaling dimensions in the 6'-plane is bounded by the parabola depicted in Fig. [2J 



Figure 2: Left figure: The points in the shaded region correspond to the scaling dimensions 
{6',6") realized by the "n_|_-family" . The bold-faced part of the vertical axis represents 
the scaling fields with n"^_ = 0. The points above the parabola satisfy n'_^ < 4. Right 
figure: The condition n'_^ > 2 is imposed in addition. It is satisfied by all points below 
the upper parabola. 

The "n_-family": There exists a scaling field (y^, T) = (0, of complex scaling 

dimensions 6 = 6' -\- iO" for every point {9' , 9") in the 6'-plane. The associated exponents 
n' + in" are given by the eqs. fl5.46p with the lower sign. 

In either family scaling fields with 9' > 0, ^' = 0, < are relevant, marginal, and 
irrelevant, respectively. 

We also found the exceptional solution fl5.44p with 9' = au"^ ^ 1.2, 9" = 0. Since 
uj ^ 1.69 it grows more slowly than ip'^ asymptotically and defines a relevant scaling field 
therefore. 

All scaling dimensions come in complex conjugate pairs. If {yg, T) is a complex 
solution of the linearized flow equation for {9', 9"), so is {y*, T*) with {9', —9"). Forming 




-fli 
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(5.50) 



linear combinations we obtain the following two real solutions: 

= e\h\ (^"' e"^'* cos(n"lnv?-^"t + 5) 

Im| ] =6\h\ 1^] e-''' sm(n"\n^-9"t + 5). 

Here we wrote T{(p) = hip"' and allowed for a complex amplitude h = \h\ e"^. In (15.501) 
the complex exponent (n', n") stands for either (n'_,_, n'^) or (n'_, n'!_) and is given in terms 
of {6', 6") by the expressions (I5.46p . Obviously all linear solutions with n" ^ are non- 
polynomial in (p, and their dependence on both ip and the scale is oscillatory. When cp 
approaches the singular point of vanishing metric, = 0, the solutions oscillate infinitely 
rapidly. Solutions with n" = do not oscillate, but are still non-polynomial unless n' G IN. 

(b) Scaling fields with t/g ^ 0. Now we switch on the inhomogeneity in eq. (15.351) . 
the terms proportional to yg. The general solution of the inhomogeneous equation is of 
the form T = Tin + '^hom where Tin is a special solution of the inhomogeneous, and Thom 
the general solution of the homogeneous equation (with the same parameter 6). Making 
a <yj^-ansatz one easily finds the following special inhomogeneous solution: 

Ti„(^) = ^+y.^3I3I^^- (5-51) 

The second condition to be satisfied is (15.341) . Inserting T = Tin + T'hom with (I5.5ip we 
are led to 



with 

Svry, (1-2A.)' 
and the ^-independent constant 



V?i ' 1 — 2 A* ip\ 



(5.53) 



K ^ A ^* (l + ^y = 2887r (1 + 2 A.) . (5.54) 
(1-2AJ^ 
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In f l5.53p it is implied that (/^i — oo, as before. The equation (15.521) determines the values 
of 9 for which scaling fields with 1/3 7^ exist. It can be solved in the following way: 

(i) We assume that its solutions 6*1^2 are such that t^hom = for the function Thom 
solving the homogeneous equation fl5.39p with 6 = 61^2- 

(ii) We determine 2 by solving the simplified equation in which "i^hom = 0: 

(iii) We prove selfconsistency by showing that the values 61^2 obtained do indeed belong 
to a homogeneous solution with ^9hom = 0. 

Eq. (15.551) is a quadratic equation for 6. Recalling the definitions of the various 
constants one finds that the solutions are a complex conjugate pair 6*1^2 = 6*' ± iO" with 

e' = A 

(5.56) 



0" = 2^2 VI + 3 -21/3 ^6.1837 
As for step (iii), the general homogeneous solution is given by {h± G C) 

Thorn = /^+^"+('''") + /^-^"-('''"^ (5.57) 

Using (I5.46P we find that for the critical exponents (15.560 the real parts of the exponents 
n+ and n_ are given by 

n+(0', e") ^ 3.87, n'_{e', 9") ^ -0.488 (5.58) 

Since both of them are smaller than 4, we indeed obtain t^hom = when we insert Thom 
into (I5.53P and let ipi 00. This proves the consistency of the procedure. 

We observe that the scaling field {yg 7^ 0, Tin) is the only one which grows oc 
asymptotically; all the others have a weaker growth and are subdominant for ip —>■ 00. 
This complex scaling field, or the two real ones equivalent to it, is exactly the one which is 
accessible to the Einstein-Hilbert truncation. In fact, the numbers (I5.56P are exactly the 
critical exponents we found in the CREH approximation [I]. They describe the spiraling 
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of the trajectories near the NGFP on the {g, A)-theory space. In the present context this 
2- dimensional space is to be regarded as a subspace of the infinite dimensional ((7, F( ■ ))- 
theory space. 

To summarize, we can describe the linear flow near the NGFP on y( ■)) -space 
as follows. There are two scaling fields with a non- vanishing (7-component. Their only 
other non- vanishing component is in the (y9^-direction of the function space {Y{ ■ )} which 
we identify with the A-direction of the CREH truncation. Those two scaling fields are 
relevant, have the critical exponents 6' = 4, 6" ^ 6.18, and coincide exactly with those of 
the CREH approximation (cf. the discussion in Subsection 14.31) . All other scaling fields 
have vanishing (7-component and correspond to perturbations of the potential alone. There 
are both relevant and irrelevant fields of this type; generically they correspond to non- 
polynomial oscillatory functions of ip. 

5.5.2 UV Critical Manifold and Subsidiary Conditions 

In the simplest case, when the eigenvectors of the stability matrix at the NGFP 
form a complete set, the subset of eigenvectors belonging to eigenvalues —6 with Re6' > 
spans the tangent space to the UV critical manifold =5^uv- Its dimensionality Auv equals 
the number of scaling fields with Re 6' > then. They are relevant in the sense that they 
grow when k is lowered. Since every complete trajectory inside ^uv defines a possible 
asymptotically safe quantum theory, there exists a Auy-parameter family of such theories 
if the solutions of the linearized flow equations generalize to solutions of the full nonlinear 
equations which extend down to /c = 0. If this latter condition is not satisfied for all 
scaling solutions the number of free parameters is smaller than Ayy. 

What the above stability analysis shows is that, to be precise, for a specific definition 
of the function space {Y{-)} there are infinitely many directions in the truncated theory 
space along which the NGFP is UV attractrive at the linearized level. Whether this result 
is directly relevant for the dimensionality of ^uv in full quantum gravity and its degree 
of predictivity is not clear yet. The actual number A of free parameters the quantum 
theory has is decided at the nonlinear level only. This number can well be different 
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from Auv which refers to (the tangent space to) J^uv at the NGFP. Even in the full 
(untruncated) theory it can happen, for instance, that not all trajectories starting there 
can be continued to = 0, whence A < Auv- Moreover, even if they reach k = 0, 
most trajectories will in general be unacceptable according to additional physical criteria 
(existence of a classical regime, etc.) which also can decrease A. This issue could be 
analyzed only by a comprehensive numerical analysis of the partial differential equation 
for Yfc. We shall not embark on this analysis here, in particular as it is not clear how 
the number A of the conformally reduced theory relates the corresponding number in the 
space of functionals depending on the full metric. 

Let us now discuss the other issue mentioned above, the precise definition of theory 
space. Already at the linear level there exists an ambiguity which has an enormous impact 
on the number of relevant directions one finds, namely the choice of the function space 
{Y{ ■ )} in which all admissible potential terms and hence scaling fields are supposed to 
"live" . In the above analysis the only condition which we imposed (besides differentiabil- 
ity) was an asymptotic growth not faster than ^9^. To show that the precise specification 
of the space {^(O) is crucial we shall now consider various plausible choices. For the 
time being we do not know which one will be ultimately correct (in the sense of closest 
to full QEG). 

(a) Subsidiary condition = 0: n' > 2 

In Subsection 14.41 we saw that the RG flow does not generate contributions to Y^^'^^ if 
^fe™^ = initially: Y^^^^ = =^ dkY^™^ = 0. Furthermore, in Subsection 15.41 we pointed 
out that, for both topologies considered, the fixed point potential has no singular part: 
y^smg _ This mcaus that the space of functions {Y \ y^™s = 0} = 3^ contains the fixed 
point and is invariant under the RG flow: No trajectory starting on it will ever leave it. 

As a consequence, we may consistently impose the subsidiary condition Y^^^^ = 
and define the theory on the smaller function space y. Imposing the condition y'^^s = q 
implies a corresponding condition on allowed basis vectors in the tangent spaces of 3^, in 
particular at the NGFP: Allowed scaling fields have vanishing components in the "sing" 
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directions. For the scaling fields of the type {ug = 0, T = this condition is met when 
T" = n{n — l) y?""^ does not blow up in the limit — > 0, i. e. when n' = Ren > 2 and in 
the special cases n = and n = 1. 

Imposing n'^ > 2 on the (0, y9")-scaling fields has the following consequences: (i) 
The ?T,_-family is eliminated completely; no G C satisfies n'_{6',6") > 2. (ii) As for the 
n+-family, n'_^_{6',6") > 2 is satisfied if 

e'<bi + b2e"^ (5.59) 

with 6i = 2(1 - 21/3/3) ^ 1.160 and 62 = 1 21/3/(2^3 - 1)2 ~ 6.217. The scaling 
dimensions {6', 9") of allowed scaling fields correspond to points in the 5-plane which lie 
below the parabola 6' = hi + 62 O""^ and, because of the n'j^ < 4 constraint, above the 
parabola fl5.49p . The corresponding region is shown in the right diagram of Fig. [2l 

(b) Integer exponents: n G !Z 

As another example, we assume that {y?", n G Z} is a "basis" on {Y{ ■ )} so that every 
function can be expanded in a Laurent series: Y{ip) = J2nez V^"- The i/g = 0-scaling 
fields are still the monomials ip"', but n is restricted to be real and integer now. As a 
consequence, the associated scaling dimension 6{n) is also real. It is given by eq. (15.401) . 
The function 6 = 6{n) for n real is plotted in Fig. [31 
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Figure 3: The function 6{n) of (15.401) for n real. It is positive only between n = and 
n = 2u^ 3.38. 



49 



(c) Real exponents: n G R 

As a last example, let us assume that for some reason the oscillatory scaling fields are to be 
excluded so that we must impose n" = 0. Using (15.461) we can show that the scaling fields 
T oc (f"" with n'J_ = have dimensions {{6', 6" = 0) \ — oo < 6' < au^} and real exponents 
n = n'^ = LU [1 ± 1 — 6' /auj^ ]. In the complex ^-plane the allowed dimensions lie on a 
part of the ^'-axis. If, on top of ri" = 0, the constraint n' < 4 is imposed we find: 

• In the ?T,_-family, all scaling fields with n"_ = also satisfy n'_ < 4. 

• In the ?T,_|_-family, the scaling fields which satisfy both = and < 4 have 
scaling dimensions 

{{e',e" = 0) I -ai <e' <auj^} 

and real exponents n = n'_^_ = lu [1 ± a/1 — Q' ja.uP' ]. In Fig. [2] these scaling fields 
correspond to the bold-faced interval on the vertical axis. 

For the time being it is not clear how to choose the function space {y( ■ )} optimally, 
but it is obvious that the count of relevant directions crucially depends on it. The optimal 
choice would be the one for which the LPA mimics full quantum gravity as well as possible. 
In order to find it one would have to gain a better understanding of the invariants of the 
full metric g^J,y^ their importance for the fiow, and their relation to the terms that could 
possibly appear in the conformally reduced action. 

We close this subsetion with several remarks. 

(i) Also the NGFP found by the method of symmetry reduction [28] where only metrics 
with two Killing vectors are quantized has infinitely many relevant directions. In fact, all 
scaling fields accessible to this approximation are relevant. Their relation to the invariants 
of the full 4-dimensional metric is likewise unclear. 

(ii) For the unconstrained {y(-)}-space the real form of the yg = 0-scaling field is 
given in (15.501) . These solutions are non-polynomial in (if ^ W). In a sense, they 
are analogous to the so-called Halpern-Huang scaling fields [62, 63] in standard scalar 
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matter field theories. Halpern and Huang analyzed their Gaussian fixed point by means 
of the LPA and found, quite unexpectedly, infinitely many relevant directions. They 
had been overlooked by perturbation theory because the eigenpotentials depend on the 
field in a non-polynomial way. As to yet, their status is still somewhat unclear. In the 

literature [64] arguments have been put forward suggesting that one should eliminate 
them by an appropriate choice of the function space. Likewise discarding them in gravity 
we are back to only three relevant directions, </?, (/'^j f^- 

(iii) The above analysis apphes to the topology. For it is not possible to obtain 
comparably transparent results in analytic form since the fixed point potential is known 
only numerically. However, we do not expect qualitatively new features to arise for S^. 

6 Transition to a Phase with Unbroken 
DifFeomorphism Invariance 

An exploration of the RG fiow and in particular the structure of =5^uv the nonlinear 
level has to rely upon numerical methods mostly. There are essentially two strategies for 
finding trajectories in c5^uv- The first one consists in a downward evolution (i.e. towards 
smaller values of k) starting at plus a relevant scaling field. We tried this method for 
R'^ where all scahng fields at the NGFP are known exphcitly; it turned out however that 
it is very difficult to implement it in a numerically reliable way because of the densely 
winding spirals near the fixed point. 

The second strategy consists in guessing an initial point in theory space and evolving 
upward (towards larger k) from this point. If the point happens to lie on ^uv the 
trajectory will hit the fixed point for k ^ oo. In practice one will vary the initial 
point until the trajectory stays close to the NGFP for a very long RG time. Small 
errors in the guessed initial point will always cause the trajectory to run away from the 
NGFP ultimately. However, by finetuning the initial point one can construct a numerical 
approximation to a trajectory in =i^uv any desired level of accuracy. By systematically 
changing the initial point it is possible to trace out the nonlinear structure of =5ijv in this 
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way, at least in principle. A complete numerical analysis of this kind is beyond the scope 
of the present paper. We shall rather present some typical trajectories found by this trial 
and error method and describe their general properties. 

Before turning to the numerical solutions we discuss a rather special, but instructive 
class of trajectories which can be found analytically. 



6.1 An Analytic Solution {R^ case) 



In the case of the topology there exists for every given trajectory of the Einstein- 
Hilbert truncation (here always taken of type Ilia) a 1-parameter family of trajectories 
in ^uv which can be found in closed form: 

1 



n(¥^) = Ci|v^ 



6 



(6.1) 



It is easy to see that (16.11) solves eq. (I3.44p . Here Ci is a free parameter with the dimension 
of a mass. It is obvious that these trajectories lie inside =5^uv- for k oo the potentials 
(16. ip approach the fixed point —A* ^p'^/Q. Consistent with the linear analysis the <^-term is 
relevant, it grows when k is lowered. (Eq. (16.11) is valid beyond the linearization though.) 
In dimensionful units the corresponding potential has a /c-independent 0-term: 



An 



6 G, 



(6.2) 



If Ci > the functions —Yk{(f) and Uk{4>) have a critical point, a global minimum, at 
nonzero values of the field: 



ipo{k) 



Mk) 



3^ ^1' 

2 ^ Afc k 

^ ^ Gk 
2^'Tk 



1/3 



1/3 



Vo{k)/k. 



(6.3) 



(6.4) 



If Ci < both —Yk{ip) and Uk{(f)), for any k, assume their global minimum at vanishing 
ip and 0, respectively. 

In the case Ci > 0, foik) vanishes only in the limit k ^ oo where foik) oc k^^^^. 
While we lower k from infinity down to the turning point scale k^, the minimum (fo{k) 
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moves towards larger yj-values. Slightly below the type Ilia trajectory enters the 
classical regime in which G^, A/^, and 0o(^) are approximately fc-independent. In this 
regime v'o(^) is heading for smaller values of (p when k is lowered further: ipo{k) = 
k4>o{k) oc k. 

Fig. m displays the exact trajectory (16. ip for Ci = +1; the underlying trajectory 
{gk, Afc) of the Einstein-Hilbert truncation was taken to be the type Ilia trajectory with 
(7t = 10~^^ at the turning point. We picked this tiny value because it leads to a long 
classical regime, see [44] and [38]. The first four plots in Fig. H] show the potential —Yk{(p) 
for different values of k. The plots (a) and (b) correspond to scales above the turning 
point, knc < k < oo, while the diagrams (c) and (d) refer to scales between the turning 
point and the termination scale, k^ > k > A;term; cf. Fig. 111(e). Here and in the following 
dimensionful quantities are always plotted in units of /ct, and the dotted (dashed) line 
represents the fixed point (turning point) potential. The gray level of the potential curves 
changes from black to gray for decreasing values of k. We always plot minus Yk{'^) 
because it is this function that corresponds to the actual potential Uk{4>) by a rescaling 
with positive factors. The plots Hl^b) and (d) are "zooms" into the small-y? region of the 
diagrams (a) and (c), respectively. We observe the following pattern which is generic for 
Ci > 0: When k is lowered from "fc = oo" the potential —Yk{f) immediately assumes its 
minimum at a nonzero value. Between k = oo and k = k^ this minimum moves towards 
larger yj- values, while it returns to smaller values once k has passed fcx- It is instructive 
to look at the latter regime in dimensionful units. Fig. Hl^f) displays the position of the 
minimum (f)o{k) = (po{k)/k below the turning point. We see that 0o is perfectly constant 
all the way down to the termination scale. (It is at about /ctorm ~ 10~^ for the trajectory 
chosen.) 

It is tempting to interpret this RG trajectory as describing a kind of continuous 
("second order") phase transition with respect to the scale, taking place for A; oo: 
At A;"="cx3 the potential —Yk{(p) has its global minimum at ^9 = 0; for all < oo it is 
located at a nonzero (p = ifo{k) ^ 0. Likewise the minimum (f)o{k) of Uk{(f)) continuously 
approaches zero for — > oo and is nonzero for all < oo. The significance of those 
minima becomes clear when we search for (running) solutions of the (running) effective 
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-J'fc(^f) decreasing k -Yk((p) decreasing k 




(f) 



Figure 4: The trajectory discussed in Subsection 16. 1[ The diagrams a and b (c and d) 
show the potential for different values of k above (below) fcx? as indicated in the schematic 
diagram e. The dotted (dashed) line is the fixed point (turning point) potential. The plot 
b (d) zooms into the diagram a (c) at small if. The plot f shows the dimensionful minimum 
00 as a function of A; < /ct- Dimensionful quantities are in units of k^- The gray level of 
the various curves changes from black to gray for decreasing values of k. 
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field equations, i.e. look for stationary points of F^. 

In fact, the vacuum expectation value of the metric can be found by solving the 
effective field equation 5Ti^/5(f) = 0. Therefore, within the LPA, the fc-dependent metric 
expectation value is {g^u)k = 0o(^)^ ^fiu- If v^o = this expectation value vanishes ac- 
cording to the LPA. Hence the phase the system is in at fc"="oo is a phase of unbroken 
diffeomorphism invariance [53,59]. In this phase the metric has a vanishing expectation 
value, with non-trivial fluctuations about {gf_iu)k = 0, but no "metric condensate" sponta- 
neously breaks diffeomorphism invariance. For k < oo a nonzero (and within the present 
approximation therefore necessarily non-degenerate) metric {g^u)k spontaneously breaks 
the diffeomorphism group down to the stability group of the ground state metric. For 
{9nv)k = <Poiky g^^ with g^^, the fiat metric on this is the (Euclidean) Poincare group. 

If 00 is approximately /c-independent over a wide range of /c-values the "frozen out" 
metric {g^i,)k = <Po defines an approximately classical spacetime. In the above example 
we saw that this is indeed the case for k between /ct and /ctorm- We also know [38, 44] 
that the classical regime predicted by the Einstein-Hilbert truncation is the longer the 
closer the turning point is to the GFP. If 0o does not depend on k we may introduce 
units such that 0o = 1 so that, with the reference metric ^f^i, = 5^^, the expectation 
value equals {g^u)k = 5^lv within the classical range of fc- values. In this way the emergence 
of an approximately classical, flat spacetime can be understood as a condensation and 
symmetry breaking phenomenon similar to those in standard scalar field theories. 

While we are familiar with gravity in the classical regime of the broken phase, it 
seems to be much more the exception than the rule. In the present setting it requires a 
very smal Rvalue of (^t and a positive constant Ci in (16.11) . For Ci < the global minimum 
of —Yk{'^) stays always at y^o = 0. So gravity is always in its "symmetric phase" and the 
metric develops no expectation value, let alone a scale independent one. 



^If one tentatively applies the discussion to the real Universe one has « 10 [40,44]. 
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-Ykitp) decreasing k -Yk(<p) decreasing k 




Figure 5: The trajectory discussed in Subsection 16. 2[ The diagrams a, b, and c show the 
potential for different values of k above k^, the diagram d for /e- values below A;t- The 
dotted (dashed) line is the fixed point (turning point) potential. The plots b and c are 
zooms into the diagram a at smaller (p values. 

6.2 First order-type Transitions {R^ case) 

In Fig. [5] we show a trajectory in ^uy, for i?^, which was obtained numerically by 
the trial and error method. The pertinent Einstein-Hilbert trajectory has = 10^^. 
The initial scale for the upward integration of the partial differential equation was chosen 
as /cinti = k'x/2 where k^ = 1 in computer units. Assuming that the initial potential ^fcinti 
is a polynomial of the form + yj^ + yj^ we finetuned the coefficients of the y?- and the 
93^-term in such a way that Yki^p) gets close to Y^{ip) and stays there for a long RG time. 
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The successful initial potential reads 

y^..^ = -0.10 V + 0.20 - i Afc,„, ^\ (6.5) 

To obtain the complete trajectory we also evolved (16. 5p downward, from A;inti to fcterm- 
The Figs. |5](a), (b), and (c) show the running potential for various values of k above k^, 
while Fig. EJ^d) shows representative potentials with k below k^. In Fig. EJ^a) we see that 
Yk{ip) approaches the (dotted) fixed point potential for large values of k. We verified that 
the approach to Y^:{ip) is oscillatory, as a consequence of the spirals in the g-X-plane. The 
Figs. [5](b) and (c) show that at a certain value of k the potential first develops an inflection 
point and then a new local minimum and maximum. At an even lower, "critical" scale 
kc the new minimum becomes the global one. Above k^ the new global minimum moves 
towards larger values, then towards smaller ones, see Fig. [n](d). The latter regime is the 
classical one again, with an approximately constant 0o = fo/k. Making smaller one 
can obtain an arbitrarily long (in k) regime with an essentially classical, flat spacetime. 

As compared to the analytic trajectory of the previous subsection two features of 
the phase transition are new here: it happens at a finite transition scale kc < oo, and it 
is discontinuous. For k > k^ the global minimum is at (po = 0, at k = k^ it jumps to a 
value 00 > 0. In analogy with the phase transitions in standard scalar theories we call 
this a first order transition. 

According to this trajectory the expectation value {g^u)k vanishes for all k > k^. 
One might wonder therefore whether the field modes with momenta above k^. are correctly 
integrated out. The answer is in the affirmative. To see why we recall that a RG trajectory 
is a running functional defined over a certain space of functions {0( • )}. The global 
minimum of is merely a single point in this space, 0o( ■ )• The flow equation has the 
structure ''kduV k[4>{ ■ )] = RHS(0(-)) for all 0( ■ )" so the point 0o(') plays no special 
role. Only after having constructed and solved the flow equation which involves arbitrary 
"off-shell" arguments 0( ■ ) we can go "on-shell" , i. e. look for solutions to the effective field 
equation ((5rfc/(50)[0o( ■ )] — 0- Only here it is decided whether the metric has a nonzero 
expectation value or not. 

In the general case of the full functional ^k[gfiu,guuy ' ' '] there exists a subtlety. 
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however, which is absent in the present setting for conformally reduced QEG. In the 
original construction of the full gravitational average action [7] it is assumed that g^^ and 
'g^^, are non- degenerate because the inverses g^'^ and g^'^ are needed there. Hence, at best, 
dfiu = and/or = are singular boundary points of the space {g^ui ■ ),g^ui ')^' ' '} 
over which Tk is defined. As a result, the original average action approach can describe 
the symmetric phase of gravity only as a singular limit. (In the full Einstein-Hilbert 
truncation, for instance, {g^u)k for k — ^ oo.) The problem of generalizing the original 
construction to a theory space \^r[g^^,g^j^, ■ ■ ■]} where T is defined for non- invert ible g^i, 
and ^^j,, too, will be addressed elsewhere. Here we only mention that in the present 
reduced setting this problem does not arise since we may use the reference metric ^^i, 
which is always non-degenerate in order to formulate the flow equation. 

A similar analysis of spontaneous symmetry breaking but in a perturbatively renor- 
malizable model of quantum gravity has been performed in ref. [53]. The scaling proper- 
ties of the potential and the average metric are quite different there as the theory is not 
asymptotically safe. In particular {g^v)k increases rather than decreases for k — > oo. 

6.3 Phase Transitions in the case 

In the previous subsection we observed that in the i?^ topology a phase transition 
can occur at finite k. It is easy to see that, in the case of S"^, a phase transition, if it 
occurs at all, must take place at a finite value of k. For 5*^ the fixed point potential 
—Y^{ip) has the shape plotted in Fig. [TJ Contrary to the simple — K oc ip^ on this 
function is structurally stable in the sense that a smooth infinitesimal deformation of — 
cannot give rise to a new local or global minimum. This was different for — oc where 
an infinitesimal term oc if can shift the global minimum from (/Jq = to y^o 7^ 0. 

As for phase transitions on S''^, both a first order and a second order scenario are 
possible. Assume, for instance, we start the downward evolution from the — pictured 
in Fig. [T] with an infinitesimal relevant perturbation added. Then it might happen that 
for a certain period of RG time, for k^ < k < 00, say, no critical points form at > 0, 
but that the slope of — Y'fc(v^) at the origin decreases from the originally positive value 
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— F/(0) > to zero at k = k^. —Yl^iQi) = 0. If the slope continues to decrease for k 
below ko then a new global minimum forms at v'o(^) > 0, and v'o(^) moves away from 
the origin in a continuous way. This, then, is a "second order" phase transition taking 
place at k = k^. The scale k^ is necessarily finite (/cq < oo) since the initial slope — Y'J(O) 
is strictly positive; hence a nonzero running time is needed in order to reduce it to zero. 
One can find trajectories which indeed realize this scenario. 

A possible first order scenario proceeds as follows. During the downward evolution 
from —Y^{ip) at A;"="oo first an inflection point forms at some y^infi > and then a new 
local minimum and a maximum arise from it. Upon further downward evolution the new 
minimum becomes the global minimum at a critical scale k^. This behavior corresponds 
to a discontinuous phase transition; the global minimum suddenly jumps from y^o = to 
7^ 0. Also in this scenario the phase transition happens at a finite scale k^ < oo. The 
reason is that the slope of —Y^ is strictly positive everywhere. Hence, again, we need a 
nonzero running time until ^^(v^) vanishes at y^infl ^ 0. 

In Fig. [6] we display a trajectory, found numerically by the trial and error method, 
which realizes this kind of first order scenario. The respective Einstein-Hilbert trajectory 
has = 10~^. At the initial scale for the upward integration of the partial differential 
equation, fcinti = 10 k^, we chose an initial potential Yk.^^^ of the type ip'^ + ip"^ + ip^ + ip^. 
After finetuning the <y9°, (p"^-, and (y9'^-coefiicients so that Yk{ip) gets close to Y^{(p) and 
stays there for a very long RG time we obtained approximately 

Yk,^M = 0-3845 + 0.9 _ o.OOl - | A,,^,^ (6.6) 

As for the smoothing scheme, the spectral functions used in the numerical integration 
were p{(p) = 1 + j^(p'^ and p{(p) = (p^. 

7 Summary and Conclusions 

In this paper we explored the "background independent" renormalization group flow 
of the effective average action for QEG on an infinite dimensional theory space. Consid- 
ering conformally reduced gravity we quantized only the fluctuations of the conformal 
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Figure 6: The 5''^ trajectory discussed in Subsection I6.3[ The diagram a shows the 
potential for different values of k above /ct- The plot b zooms into diagram a for small 
if and shows the first order phase transition. The dotted (dashed gray) line is the fixed 
point (initial) potential. The gray level of the various curves changes from black to gray 
for decreasing values of k. 

factor, employing the Local Potential Approximation for its effective average action. The 
relative simplicity of the resulting system of flow equations, a partial differential equa- 
tion coupled to an ordinary one, allowed us detailed investigations which at present are 
prohibitively complicated in the full theory. While at first sight reminiscent of the RG 
equations for standard scalar theories, the requirement of a "background independent" 
quantization results in crucial differences. In the infinite dimensional space of potential 
functions we found both a Gaussian and a non-Gaussian fixed point; they generalize the 
fixed points known from the (conformally reduced) Einstein-Hilbert truncation. The scal- 
ing dimensions at the GFP were found to be very different from those in a scalar theory, as 
a consequence of the built-in "background independence" . The results on the structure of 
the RG flow provide further non-trivial evidence for the viability of the asymptotic safety 
scenario. 

We studied in detail the linearized RG flow near the NGFP and we found that for 
some choices of the space of scaling fields infinitely many relevant directions can arise. 
The corresponding scaling fields are non-polynomial functions similar to the Halpern- 
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Huang directions at the GFP of conventional scalar field theory. As to yet, their status 
and implications for the actual predictivity of the theory is not understood, not even 
in the case of standard scalars. At least in gravity the LPA is not sufficient in order 
to meaningfully address this question; in particular it is unclear whether the Halpern- 
Huang-like scaling fields can descend from invariants built from the full metric which 
occur in the fixed point action. 

In QEG the beta functions, by construction, are independent of any specific metric. 
They do depend on the topology of spacetime, however. We illustrated this point by 
considering the cases -S"^ and i?^ in parallel. 

Analyzing the general properties of the RG flow and in particular of the UV critical 
manifold associated with the NGFP we saw that on J^uv the potential Yk{ip) behaves 
as for ip — > oo. Recalling that for every term in rfc[(7^,^] which has the structure 
f d^Xy^ (curvature)" contributes to the (/^'^^^"-monomial wc can conclude that, in the 
conformally reduced theory, we see no "shadows" of invariants with negative powers of 
the curvature. Leaving aside the possibility of cancellations among invariants of the 
same dimension, this result suggests that terms of the form J d^Xy^i?"", n — 1,2,3, ■ ■ ■ , 
should not occur in RG trajectories on S^uy. In fact, recently Machado and Saueressig [21] 
have analyzed 3-parameter truncations where a single term of this type was added to the 
Einstein-Hilbert truncation. And indeed for moderate n it turned out that on this 3- 
dimensional theory space no satisfactory UV fixed point could be found. In the light of 
the present results the natural explanation which suggests itself is that the trajectories 
of [21] arc not on S^uy when the non-local invariant is included. The same remark applies 
to J d'^Xy/g' In R which produces a term (/?^ In (p which likewise grows faster than (p'^ and, 
exactly as expected, spoils the NGFP. (The actual motivation in [21] for including these 
invariants was the hope that they would improve the IR behavior, and this indeed turned 
out to be the case.) 

In the opposite hmit ^ 0, the NGFP potential Y^:{(p) was found to be regular. We 
also saw that after appropriately constraining the scaling fields even the running potential 
Yk{ip) is regular along all trajectories in =^uv- We can interpret this regularity by saying 
that Kf. contains no "shadows" of the invariants /„ = J d'^'Xy/g' (curvature)" oc (^^~^" for 
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n = 3, 4, ■ ■ ■ . The monomials Iq o: (p^ , Ii (x (p^ and I2 = const are regular; they represent 
the cosmological constant, the Einstein-Hilbert term, and the (curvature) ^-terms, and 
they have a nonzero fixed point coupling^. For the higher ones I3 cx 9?^^, /4 cx v?"^, ■ ■ ■ 
the fixed point coupling is zero. 

Here we see that in the conformally reduced framework the Goroff-Sagnotti term 
/ d'^x^/g' R^uap R'^^pa R^'^^^ which contributes to y?"^ does not play any special role for 
the UV behavior of the theory. Imposing the subsidiary condition Y^^'^^ = it is even 
absent all along the trajectories in =5^uv- This term is known to spoil the perturbative 
renormalizability of full quantum gravity at the two-loop level [58]. According to the 
general theory [29] and well known examples [65] perturbative nonrenormalizability can 
well coexist with asymptotic safety. The present results are fully consistent with this 
picture. 

Solving the fiow equation for the potential numerically we obtained examples of RG 
trajectories in the NGFP's ultraviolet critical hypersurface ^uv- The quantum theories 
based upon some of them show phase transitions from a phase in which the global mini- 
mum of Uk{(f)) is at some 7^ to a phase where it is at = 0. In the former situation 
the conformal factor has a nonzero expectation value, while it vanishes in the latter. We 
interpret this as an indication that at some large (but in some cases finite) value of k 
there occurs a transition from the familiar low energy phase of gravity with a nonzero 
expectation value of the metric to a new phase in which g^^^, = (7^^) = 0. In the "low 
energy" phase diffeomorphism invariance is spontaneously broken down to the stability 
group of the metric condensate (7^;^) 7^ 0, while it is unbroken in the "high energy" phase. 
There also exist trajectories along which the symmetry is always unbroken, i. e. the met- 
ric never develops an expectation value. The corresponding quantum theories have no 
classical regime. 



^The /2-coupling , the constant piece in Y"*, is nonzero at least for S''*, see Fig. [TJ 
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APPENDIX 



A Generalizations for d Spacetime Dimensions 

In this appendix we tabulate the key formulas of the main text generalized for an 
arbitrary spacetime dimensionality d. 

The LPA truncation ansatz generalizing the average action (13.41) is given by 



2tt {d - 2) Gk 



(A.l) 



+ f-^Ffc(fV- (XB + /))} 



where u = 2/{d — 2). (Cf. ref. [22].) The running potential of the CREH truncation reads 



Fr^"(x) = CO 



d{d-2) 



X 



(d-2)p 



d-2 



The cutoff operator IZk must be chosen as 

d-l 1 
"27r {d-2) Gk 



SO that the cutoff action (13. lip becomes 

d-l 1 
^An {d-2)Ch 



AkS[f;xB,gf,u] 



In terms of the new field variables 



(A.2) 



(A.3) 



Xl^e Id'xVd f{x)R^'\-^)f{x). (A.4) 



Xb = r^'"" Xb, 



(A.5) 



the generating equation for the RG flow on 5'°' reads 

Hid i7W(-e)7M.«.^^(^ 



Tr 



i?(o) 



□ 



Xb 



1.2 

R ft/ 



□ 



□ 



-□ + Xb 



□ 



Xb 



21/ 1.2 



+ F;!{xB + f{x)) 



(A.6) 
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Eq. (1A.6P reduces to fl3.2ip in 4 dimensions. 

The negative eigenvalues and degeneracies of the Laplace-Beltrami operator on a 
unit d-sphere [60] are, respectively, 

, ^ (2n + d-l)(n + d-2)\ , , 

£n;d = n{n + d-l) and D^-d = ^, (A.7) 

The spectral functions of fl3.30p generalize to Pd{^) = Yl'^=o ^n;d 0{v'^'^~£n;d) = -'d(^max(v^; d)) 

and Pd(v9) = Y.n=Q^n;dDn;dd{v'^'' " £n;d) = <^d(^max(V5; d)) with the SUUIS 



^ {2n + d~ l){n + d-2)\ 



^ ' n=0 



^ 'n + rf-l)(2n + rf-l)(n + (i-2)! 



(A.8) 



n=0 

where = rijaaJy^P', d) is the largest positive integer satisfying S^-d = N{N + d— 1) < ip'^'^. 
For (f ^ 1 we may approximate Pd(v5) ~ ^ v^'^'' and pd(v^) ~ "^(dt^)] '{'^'^^'^^^ which boils 
down to (13.351) for 4 dimensions. 

The partial differential equations for the potential Fk{4>) and its dimensionless analog 
Yk(Lp) = k'^"'^ Fk^ip/k^^^) read in arbitrary dimensions, respectively, 

[d - 1) ad k^cj)^" + F^W) 

and, with (/9 = k^^'^ 0, 

27r(rf - 2) {1 - f } + f ^^-^^^ 



Their counterparts in 4 dimensions are given by eqs. (13.281) and fl3.4ip . respectively. Here, 
(Trf = 27r('^+i)/Vr((rf+l)/2). 

With the optimized shape function the anomalous dimension r/N for the topology 
assumes the form 
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which reduces to fl3.52p in 4 dimensions. In terms of the dimensionless quantities Qk and 
Ffc the equation flA.lOl) reads, with Lpi = k^^" cpi, 



m{9k, [Yk]) 



^71 Vd {d-2) 



(d+2)u 



d{d-l) ^'-^^ [^2. + 

with the usual abbreviation Vd = \2'^^^ tt"'/^ r((i/2)] \ 



(A.ll) 



65 



References 

[1] For a general introduction see C. Kiefer, Quantum Gravity, Second Edition, 
Oxford Science Publications, Oxford (2007). 

[2] L. Rosenfeld, Ann. der Physik 5 (1930) 113; Z. fiir Physik 65 (1930) 589; 
M. Bronstein, Phys. Zeitsclirift der Sowjetunion 9 (1936) 140. 

[3] A. Aslitekar, Lectures on non-perturbative canonical gravity, 
World Scientific, Singapore (1991); 

A. Ashtekar and J. Lewandowski, Class. Quant. Grav. 21 (2004) R53. 

[4] C. Rovelli, Quantum Gravity, Cambridge University Press, Cambridge (2004). 

[5] Th. Thiemann, Modern Ganonical Quantum General Relativity, 
Cambridge University Press, Cambridge (2007). 

[6] S. Weinberg in General Relativity, an Einstein Gentenary Survey, 

S.W. Hawking and W. Israel (Eds.), Cambridge University Press (1979); 
S. Weinberg, |hep-th/9702027] 



[7] M. Renter, Phys. Rev. D 57 (1998) 971 and hep-th/9605030 



[8] D. Don and R. Percacci, Class. Quant. Grav. 15 (1998) 3449. 



[9] O. Lauscher and M. Renter, Phys. Rev. D 65 (2002) 025013 and hep-th/0108040 



[10] M. Renter and F. Saueressig, Phys. Rev. D 65 (2002) 065016 and hep-th/01 10054 



[11] O. Lauscher and M. Renter, Phys. Rev. D 66 (2002) 025026 and hep-th/0205062 



[12] O. Lauscher and M. Renter, Class. Quant. Grav. 19 (2002) 483 and hep-th/01 10021 



[13] O. Lauscher and M. Renter, Int. J. Mod. Phys. A 17 (2002) 993 and |hep-th/0112089 
[14] W. Souma, Prog. Theor. Phys. 102 (1999) 181. 



66 



[15] M. Reuter and F. Saueressig, Phys. Rev. D 66 (2002) 125001 and |hep-th/0206145 
Fortschr. Phys. 52 (2004) 650 and|hep-th/0311056[ 



[16] A. Bonanno and M. Reuter, JHEP 02 (2005) 035 and hep-th/0410191 



[17] For reviews on QEG see: M. Reuter and F. Saueressig, l arXiv:0708. 13171 [hep-th]; 
O. Lauscher and M. Reuter in Quantum Gravity, B. Fauser, 



J. Tolksdorf and E. Zeidler (Eds.), Birkhauser, Basel (2007) and |hep-th/0511260 

0. Lauscher and M. Reuter in Approaches to Fundamental Physics, 

1. -O. Stamatescu and E. Seiler (Eds.), Springer, Berhn (2007). 

[18] R. Percacci and D. Perini, Phys. Rev. D 67 (2003) 081503; 

Phys. Rev. D 68 (2003) 044018; Class. Quant. Grav. 21 (2004) 5035. 

[19] A. Codello and R. Percacci, Phys. Rev. Lett. 97 (2006) 221301; 

A. Codello, R. Percacci and C. Rahmede, Int. J. Mod. Phys. A23 (2008) 143. 

[20] D. Litim, Phys. Rev. Lett. 92 (2004) 201301; AIP Conf. Proc. 841 (2006) 322; 
P. Fischer and D. Litim, Phys. Lett. B 638 (2006) 497; 
AIP Conf. Proc. 861 (2006) 336. 

[21] P. Machado and F. Saueressig, Phys. Rev. D 77 (2008) 124045. 

[22] M. Reuter and H. Weyer, Phys. Rev. D 79 (2009) 105005 and larXiv:0801. 32871 [hep- 
th]. 



[23] O. Lauscher and M. Reuter, JHEP 10 (2005) 050 and hep-th/0508202 



[24] M. Reuter and J.-M. Schwindt, JHEP 01 (2006) 070 and hep-th/ 05 11 021 



[25] M. Reuter and J.-M. Schwindt, JHEP 01 (2007) 049 and |hep-th/0611294 
[26] J.-E. Daum and M. Reuter, preprint Eiv: 0806. 3907) [hep-th]. 
[27] F. Neugebohrn, larXiv:0704.3"205] [hep-th] 



67 



[28] P. Forgacs and M. Niedermaier, |hep-th/0207028 



M. Niedermaier, JHEP 12 (2002) 066; Nucl. Phys. B 673 (2003) 131; 
Class. Quant. Grav. 24 (2007) R171. 

[29] For detailed reviews of asymptotic safety in gravity see: 

M. Niedermaier and M. Renter, Living Reviews in Relativity 9 (2006) 5; 
R. Percacci, larXiv:0709.3"85T] [hep-th]. 

[30] L.F. Abbott, Nucl. Phys. B 185 (1981) 189; 

B. S. DeWitt, Phys. Rev. 162 (1967) 1195; 

M.T. Grisaru, P. van Nieuwenhuizen and C.C. Wu, Phys. Rev. D 12 (1975) 3203; 
D.M. Capper, J.J. Dulwich and M. Ramon Medrano, Nucl. Phys. B 254 (1985) 737; 
S.L. Adler, Rev. Mod. Phys. 54 (1982) 729. 

[31] C. Wetterich, Phys. Lett. B 301 (1993) 90. 

[32] M. Renter and C. Wetterich, 

Nucl. Phys. B 417 (1994) 181, Nucl. Phys. B 427 (1994) 291, 

Nucl. Phys. B 391 (1993) 147, Nucl. Phys. B 408 (1993) 91; 

M. Renter, Phys. Rev. D 53 (1996) 4430, Mod. Phys. Lett. A 12 (1997) 2777. 

[33] J. Berges, N. Tetradis and C. Wetterich, Phys. Rep. 363 (2002) 223; 

C. Wetterich, Int. J. Mod. Phys. A 16 (2001) 1951. 

[34] For reviews of the effective average action in Yang-Mills theory see: 



M. Renter, |hep-th/9602012[ J. Pawlowski, Ann. Phys. 322 (2007) 2831; 



H. Gies, |hep-ph/061lT46l 



[35] A. Bonanno and M. Renter, Phys. Rev. D 62 (2000) 043008 and |hep-th/0002196; 



Phys. Rev. D 73 (2006) 083005 and |hep^th/0602159t 



Phys. Rev. D 60 (1999) 084011 and gr-qc/9811026 



[36] M. Renter and E. Tuiran, ,hep-th/0612037l 



[37] A. Bonanno and M. Renter, Phys. Rev. D 65 (2002) 043508 and hep-th/0106133 



68 



[38] M. Reuter and F. Saueressig, JCAP 09 (2005) 012 and|hep-th/0507167 



[39] A. Bonanno and M. Reuter, Phys. Lett. B 527 (2002) 9 and astro-ph/0106468 



Int. J. Mod. Phys. D 13 (2004) 107 and |astro-ph/02 10472 
E. Bentivegna, A. Bonanno and M. Reuter, 
JCAP 01 (2004) 001 and |astro-ph/03 03150[ 



[40] A. Bonanno and M. Reuter, JCAP 08 (2007) 024 and larXiv: 0706.01741 [hep-th]. 

[41] A. Bonanno, G. Esposito and C. Rubano, Gen. Rel. Grav. 35 (2003) 1899; 
Class. Quant. Grav. 21 (2004) 5005; 
A. Bonanno, G. Esposito, C. Rubano and P. Scudellaro, 
Class. Quant. Grav. 23 (2006) 3103 and 24 (2007) 1443. 



[42] M. Reuter and H. Weyer, Phys. Rev. D 69 (2004) 104022 and hep-th/0311196 



[43] M. Reuter and H. Weyer, Phys. Rev. D 70 (2004) 124028 and hep-th/0410117 



[44] M. Reuter and H. Weyer, JCAP 12 (2004) 001 and hep-th/0410119 



[45] F. Girelli, S. Liberati, R. Percacci and C. Rahmede, 
Class. Quant. Grav. 24 (2007) 3995. 

[46] D. Litim and T. Plehn, Phys. Rev. Lett. 100 (2008) 131301. 

[47] J. Moffat, JCAP 05 (2005) 2003; 

J.R. Brownstein and J. Moffat, Astrophys. J. 636 (2006) 721; 
Mon. Not. Roy. Astron. Soc. 367 (2006) 527. 

[48] J. Ambj0rn, J. Jurkiewicz and R. Loll, Phys. Rev. Lett. 93 (2004) 131301. 

[49] J. Ambj0rn, J. Jurkiewicz and R. Loll, Phys. Lett. B 607 (2005) 205. 

[50] J. Ambj0rn, J. Jurkiewicz and R. Loll, Phys. Rev. Lett. 95 (2005) 171301; 
Phys. Rev. D 72 (2005) 064014; Contemp. Phys. 47 (2006) 103. 



69 



[51] A.M. Polyakov, Yad. Fiz. 64 (2001) 594 

[English Translation: Phys. Atom. Nucl. 64 (2001) 540]. 

[52] R. Jackiw, C. Niiiiez and S.-Y. Pi, Phys. Lett. A 347 (2005) 47. 



[53 

[54; 
[55 
[56 

[57 



[58; 

[59 
[60 
[61 

[62 

[63 

[64 



R. Floreanini and R. Percacci, Nucl. Phys. B 436 (1995) 141; 
Phys. Rev. D 46 (1992) 1566. 

K. Symanzik, Nuovo Cim. Lett. 6 (1973) 77. 

For a historic account see: G. 't Hooft, Nucl. Phys. B 254 (1985) 11. 
O. Lauscher, M. Renter and C. Wetterich, 



Phys. Rev. D 62 (2000) 125021 and hep-th/0006099 



J.V. Narlikar and T. Padmanabhan, 

Gravity, Gauge Theories and Quantum Cosmology, D. Reidel, Dordrecht (1986), 
Chapter 12 and references therein. 

M.H. Goroff and A. Sagnotti, Phys. Lett. B 160 (1985) 81; 
A.E.M. van de Ven, Nucl. Phys. B 378 (1992) 309. 

E. Witten, Nucl. Phys. B 311 (1988) 46. 

M.A. Rubin and C.R. Ordonez, J. Math. Phys. 25 (1984) 2888; 26 (1985) 65. 

D. Litim, Phys. Lett. B 486 (2000) 92; Phys. Rev. D 64 (2001) 105007; 
Int. J. Mod. Phys. A 16 (2001) 2081. 

K. Halpern and K. Huang, Phys. Rev. Lett. 74 (1995) 3526; 
Phys. Rev D 53 (1996) 3252. 

A. Bonanno, Phys. Rev. D 62 (2000) 027701; 
H. Gies, Phys. Rev. D 63 (2001) 065011. 

T.R. Morris, Phys. Rev. Lett. 77 (1996) 1658; Phys. Lett. B 334 (1994) 355. 



70 



[65] G. Parisi, Nucl. Phys, B 100 (1975) 368; Nucl. Phys, B 254 (1985) 58; 
K. Gawedzki and A. Kupiainen, Nucl. Phys. B 262 (1985) 33; 
Phys. Rev. Lett. 54 (1985) 2191; Phys. Rev. Lett. 55 (1985) 363; 

B. Rosenstein, B.J. Warr and S.H. Park, Phys. Rept. 205 (1991) 59; 

C. de Calan, P.A. Faria da Veiga, L. Magnen and R. Seneor, 
Phys. Rev. Lett. 66 (1991) 3233. 



71 



